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INTRODUCTION 


Many techniques for solving equations involve iteration, otherwise known as the 
method of ‘refining guesses’. For example, one way to solve a real polynomial 
equation p(x) = 0 is to use the Newton-Raphson method, based on the 
recurrence relation: 


x n+1 = x r 


P( x «) 
P'{ x nY 


n = 0,1,2,... . 


The justification of this recurrence relation is shown in Figure 0.1, which 
indicates that if x n is close to a zero a of p, then x n -|_i is (usually) even closer. 
In particular, if p(x) = x 2 — 2, then the Newton-Raphson recurrence relation is 



In this unit we study the behaviour of complex iterative processes of the form 


z n +i=f(z n ), n = 0,1,2,..., 


where / is an analytic function. Usually / will be a polynomial function, 
although we spend a short time in Section 1 discussing the Newton-Raphson 
method for complex polynomials, which involves the iteration of a rational 
function. Most of Section 1 is devoted to the basic notions associated with 
iteration: nth iterates, fixed points and conjugate iteration sequences. 

Section 2 takes up the iteration of complex quadratic functions, which is the 
main topic of the unit. We begin by showing that it is sufficient to consider 
quadratic functions belonging to the basic family {P c (z) = z 2 + c: c G C}, 
because each quadratic function is equivalent, in a certain sense, to one of this 
form. We then consider, for any given c 6 C, the set of points E c which ‘escape 
to oo’ under iteration of P c , and establish some basic properties of E c and its 
complement K c . To determine some of the points of K c , we introduce the idea 
of a periodic cycle of points, and we show that certain periodic cycles must lie 
in the interior of K c , whereas others lie on the boundary of K c . This boundary 
is called the Julia set J c of P c (see Figure 0.2). 

In Section 3, we use a technique called graphical iteration, which applies only to 
real functions, to obtain various properties of the set K c when c is real. For 
example, we determine the nature of the set K c f~l IR for all real numbers c. 

In Section 4, we define the Mandelbrot set, which is the set M of numbers c 
such that the set K c is connected, and we obtain a criterion for c to belong 
to M. This criterion makes it possible to plot pictures of M (see Figure 0.3) 
and thus reveal its immensely complicated structure. We investigate this 
structure by using a result which states that if P c has a so-called attracting 
cycle, then c £ M. 


You may be tempted to think that the Mandelbrot set is in some sense an 
oddity, arising from some special property of quadratic functions. However, in 
Section 5, we describe briefly some different iteration sequences, and we find 
that the Mandelbrot set is a truly universal object. 



An initial guess of xo = 2 yields 
the sequence of refined guesses 
2,1.5,1.417,..., which 
converges rapidly to the known 
solution \/2. 


Figure 0.2 A Julia set 



Figure 0.3 
The Mandelbrot set 
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Study guide 

You may find some of the ideas in this unit quite challenging. However, your 
labours will be rewarded by gaining some insight into a fascinating subject, 
which is the object of current research. Some of this research is mentioned in 
the discussion on the audio tape, which has no associated frames and which 
may be listened to at any time. 

Section 5 and Subsection 2.4 are intended for reading only. 

Associated with this unit is a segment of the Video Tape for the course. 
Although this unit text is self-contained, access to the video tape will enhance 
your understanding. Suitable points at which to view the video tape axe 
indicated by a symbol placed in the margin. 


1 ITERATION OF ANALYTIC FUNCTIONS 


After working through this section, you should be able to: 

(a) calculate several terms of a given iteration sequence ; 

(b) determine the nth iterate of a given analytic function, for small values of n; 

(c) determine the fixed points of certain analytic functions and find their 
nature; 

(d) calculate conjugate iteration sequences', 

(e) determine the Newton-Raphson function N for a given polynomial 
function p, and describe the iterative behaviour of N when p is a quadratic 
function. 


1.1 What is an iteration sequence? 

Any sequence {z n } defined by a recurrence relation of the form 
z n +i=f{z n ), n = 0 ,l, 2 ,..., 

where / is a function, is called an iteration sequence with initial term z 0 . 
For example, if f(z) = 2z, then the recurrence relation is 
z n +i = 2 z n , n = 0 ,l, 2 ,.... 

If z 0 = i, then the first few terms of the corresponding iteration sequence are 
z 0 = i, Z\ = 2 i, z 2 = 4 i, 23 = 8 *,... . 

In this case the sequence {z n } tends to infinity, but if we had chosen z 0 = 0, 
then the corresponding iteration sequence would be constant: 

zq = 0, z\ = 0, z 2 = 0, Z3 = 0 ,... . 

Thus the behaviour of a given iteration sequence depends not only on the 
function /, but also on the choice of initial term zo■ We often represent such 
iteration sequences, as in Figure 1.1, by plotting the points zo, z \,..., and 
indicating how these points are related by the function f. Note that, in effect, 
we are plotting the domain and codomain of / on the same diagram! 


Some texts call zo the seed and 
{z n } the orbit of zo- Also note 
that the study of iteration 
sequences is often called 
dynamical systems. 



Figure 1.1 




Problem 1.1___ 

Calculate and plot the terms up to z 3 of each of the following iteration 
sequences. Also, write down the corresponding functions /. 

(a) z n+1 = z\, zo = i (b) z n+1 = \z n + 1, z 0 = 0 
(c) z n+1 = zl - 1, z 0 = 0 (d) z n+l =zl + i, z o = 0 


In order to study iteration sequences systematically, it is useful to introduce 
various basic notations. For example, if the sequence {z n } is defined by 

z n +i=f(z n ), n = 0,1,2,..., 

then 

21 = /(*o), 22 = f(z 1 ) = f(f(z 0 )), Z 3 = f(z 2 ) = /(/(/( 20 ))), 

and, in general, 

*» = /(/(•• •(/(*)))•••)), for n = 1,2,..., (1.1) 

where the function / is applied n times. We introduce a notation for such 
repeated compositions. 


For example, 
f\z) = f( Z ) 
and 

f{z) = /(/(*)). 


Remarks 

1 There is, of course, a possible confusion between 

f = f of and f = / x /, 
but the context will make the intended meaning clear. 

2 Equation (1.1), relating the general term z n to the initial term z 0 , can now 
be written in the convenient form 

z n = f n (z 0 ), n = 1,2,...; 
see Figure 1.2. 

3 Note that if m,n > 0, then 

f m (r(z)) = r +n (z) = r(r(z)), 

since composition of functions is associative. 

Example 1.1 

Determine the rules for the functions / 2 and / 3 when f(z) = z 2 — 1 

Solution 

By the definition, 

f\z) = f(f(z)) 

= f(* 2 ~ 1) 

= (2 2 -l) 2 -l 

= * 4 - 2 * 2 , 

and so 

f(z) = f(f 2 (z)) 

= {z 4 — 2 2 2 ) 2 — 1 

= z 8 - 4z 6 + 4 z 4 - 1. ■ 

6 


Alternatively, 

f\z) = / 2 (/( 2 )) 

= (* 2 -l) 4 -2(2 2 -l) 2 

= z 8 -4z 6 + 4z 4 - 1. 


There should be no confusion 
with the nth derivative f <n ^. 


!(*>)• t 

o) 

Figure 1.2 


/ 2 (z o)*'- 


Definition The nth iterate of a function / is the function obtained by 
applying the function / exactly n times: 

r = fof 0 ...of. 

Also, f° denotes the identity function f°(z) = z. 





Problem 1.2- 

Determine the rules for the functions f 2 and f 3 when f(z) = \z + 1. 


The solution to Example 1.1 suggests that for some functions / it may be 
difficult to find a general formula for the nth iterate, However, there are a 
few simple cases which are very useful. 

Example 1.2 

Find a formula for the nth iterate f n of each of the following functions. 

(a) f(z) = az, where a € C (b) f(z) = z 2 

Solution 

(a) By definition, 

f\z) = f(z) = az, 

f 2 {z) = /(/(*)) = f(az) = a(az) = a 2 z, 
f 3 (z) = f(f 2 (z)) = f(a 2 z) = a{a 2 z) = a 3 z, 
and, in general, 

f n (z) = a n z, for n = 1,2,... . 

(b) By definition, 

/ 1 (z) = /(z) = z 2 , 

f 2 (z) = f(f(z)) = f(z 2 ) = (z 2 ) 2 =z\ 

/ 3 (z) = /(/ 2 (z))=/(z 4 ) = (z 4 ) 2 = z 8 , 
and, in general, 

/ n (z) = z 2 ", for n = 1,2,... . ■ 

The following problem deals with other nth iterates which can be found 
explicitly. 

Problem 1.3_ 

Find a formula for the nth iterate f n of each of the following functions. 

(a) f(z) = z + b, where b € C 

(b) f(z) = z 3 


The formulas obtained in Example 1.2 and Problem 1.3 can be used to 
determine the behaviour of the corresponding iteration sequences, as we now 
illustrate. 


This can be proved by 
Mathematical Induction. 





Example 1.3 

(a) Prove that if f(z) = az, where |a| < 1, and z 0 G C, then 

z n = f n (z 0 ) -^Oasn-»oo. 

(b) Prove that if f(z) = z 2 and \z 0 \ < 1, then 

z n - / n (z 0 ) -* 0 as n -nx). 

(c) Prove that if f(z) = z + b, where 6^0, and z 0 G C, then 

z n = f n (zo) —> oo as n —► oo. 

Solution 

(a) By Example 1.2(a), 

z n = f n ( z o) = a n zo, for n = 1,2,... . 

Since |a| < 1, we deduce that {z n } is a null sequence, as required. 

(b) By Example 1.2(b), 

z n = f n {z 0 ) = zf, for n = 1,2,... . 

Now we use the facts that |z 0 | < 1 and 2 n > n, for n = 1,2,..., to deduce 
that 

\z n \ = \zo\ 2 < |z 0 | n , for n = 1,2,... . 

Since |zo| < 1, {|^o| n } is a null sequence and so, therefore, is {z n }, by the 
Squeeze Rule. 

(c) By Problem 1.3(a), 

z n = f n (zo) = z 0 + nb, for n = 1,2,... . 

Thus 

J_ _ 1 

z n z 0 + nb 
1/n 

z 0 /n + b 
0 

—► - = 0asn->oo (since b ^ 0); 
so, by the Reciprocal Rule, 

z n —> oo as n —» oo. ■ 

Remark The sequence {zq”} in part (b) tends to 0 extremely quickly if 
|zo| < 1. In fact, for any a such that 0 < |a| < 1, we have 



Thus {*g n } tends to 0 faster than {a n } for any a with 0 < |a| < 1. 

In Example 1.3(a) we saw that f n (z 0 ) —> 0 as n —> oo for all choices of initial 
term z 0 , whereas in Example 1.3(b), f n (z 0 ) -> 0 as n -* oo, whenever |z 0 | < 1. 
In the next problem we ask you to investigate what happens for this latter 
function / with other initial values. 

Problem 1.4___ 

With f(z) = z 2 , determine the behaviour of the iteration sequence 
z n = f n (zo), n = 1 , 2 ,..., 

when 

(a) z 0 = 1; (b) z 0 = ~i\ (c) z 0 = e 2 ^/ 3 ; (d) |z 0 | > 1. 


Unit A3, Theorems 1.2 and 1.3 


By the Binomial Theorem, 
2" = (1 + l) n 

= 1 + n + • • • > n. 
Unit A3, Theorem 1.1 


Unit A3, Theorem 1.5 


To prove this, note that if 

a n = M 2 7M n , n = 0,1,2,. 
then 



Hence {a n } is null, by Unit B3, 
Theorems 1.9 and 1.7, for 
example. 




1.2 Fixed points 

Whenever an iteration sequence, defined by a continuous function /, converges 
to a limit a, say, then the point a has the property that /(a) = a, as we now 
show. Suppose that the iteration sequence z n = f n (zo) is such that 
z n —* a as n —► oo. 

Then z n + 1 —► a as n —► oo, so that 

a = Jim^n+i = \im° f(z n ) = /(a), 

because the function / is continuous at a. As the limit a satisfies f(a) = a, it 
is called a fixed point of the function f (see Figure 1.3). 


Definition A point a is a fixed point of a function / if /(a) = a. 


For example, the function f{z) = 2z has 0 as a fixed point, whereas the function 
f(z) = z 2 has 0 and 1 as fixed points. In general, we find the fixed points (if 
any) of a given function / by solving the fixed point equation f(z) = z. 

Problem 1.5- 

Determine the fixed points of each of the following functions /. 

(a) f(z) = \z+ 1 (b) f(z) = z 2 — 2 (c) f(z) = z 3 


The behaviour of an iteration sequence z n = f n (zo), n — 1,2,..., near a fixed 
point a of an analytic function / depends to a very great extent on the 
derivative of / at a. If |/'(a)| < 1, then, to a good approximation, / maps 
small discs with centre a to even smaller discs with centre a. Thus an initial 
term zo near a gives rise to an iteration sequence which is attracted to a (see 
Figure 1.4). 


Theorem 1.1 Let a be a fixed point of an analytic function / and 
suppose that |/'(q)| < 1. Then there exists r > 0 such that 

Jirn^ f n (z 0 ) = a, for | z 0 - a| < r. 


Proof We first choose a real number a such that |/'(a)| < a < 1. Since 
f'(a) = lim and a - |/ , (a)| > 0, 

2—>cr z — a 

there is a positive number r such that 

| — z r _ - ^ a) - /'(a) | < a - imi, for 0 < \z - a| < r; 
see Figure 1.5. 

Hence 


z - a 

and so, since /(a) = a, 

l/(z)-<*| < a|z — a|, 
Thus if |zo - a | < r, then 
l/(2o) -a| < a|z 0 -a|, 
so that |/(z 0 ) - a| < r, also. 


< a, 


for 0 < \z — a| < r, 


for \z — a| < r. 


The sequence {z n+ i} is just the 
sequence {z„} with its first term 
removed. 


d a / 

2 0 * 


Figure 1.3 


Recall from Unit A4, 
Subsection 1.5, that f'(a) acts 
as a complex scale factor at a. 



Figure 1.4 


Apply the e-6 definition of limit 
(Unit A3, Section 3) with 
£ = a-|/'(a)| 
and put r = 6. 


m - /(a) 



Figure 1.5 






Hence 


\f 2 (z 0 ) — a\ < a\f(zo) - a\ < a 2 \zo - q|, 
and, in general, 

l/"(*o) - a| < o n |z 0 - a\, for n = 1,2,... . (1.2) 

Since 0 < a < 1, the sequence {a n } is a null sequence. Thus, for \z 0 — a| < r, 
the sequence {/ n (z 0 ) — a} is null, by the Squeeze Rule, and so 
lim° f n (z 0 ) = a. ■ 


In this proof, notice that the smaller |/'(a)| is, the smaller we can choose the 
number a (such that |/'(a)| < a < 1) and so, by Inequality (1.2), the faster the 
sequence {/"(zo)} converges to a. This convergence is very fast if /'(a) = 0. 

If a is a fixed point of / for which |/'(a)| > 1, then we should expect initial 
terms zq, near a (but not at a) to be pushed away from a by /. If \f'(a)\ = 1, 
then the behaviour depends in a more subtle way on the precise value of f'(ct). 
These observations suggest the following classification of fixed points. 


Definitions The fixed point a of an analytic function / is 

(a) attracting, if |/'(a)| < 1; 

(b) repelling, if |/'(a)| > 1; 

(c) indifferent, if |/'(a)| = 1; 

(d) super-attracting, if f'(a) = 0. 


For example, the function f{z) = az, where a e €, has 0 as a fixed point and, 
since f'(z) — a, this fixed point is attracting if |a| < 1, repelling if |a| > 1 and 
indifferent if |a| = 1. 

Problem 1.6____ 

For each of the following functions /, classify the given fixed point a. 

(a ) f(z) = z 2 , a = 0,1 

(b) /(z)=i z + l, a = 2 

(c) f(z) = z 2 -2, a = 2 


For any given function / with an attracting fixed point a, it is natural to ask 
exactly which points z are attracted to a under iteration of / (that is, 
f n (z) —* a-as n-*oo) and so we make the following definition. 


Definition If a is an attracting fixed point of an analytic function /, 
then the basin of attraction of a under / is the set 

{z : f n {z) —► a as n —► oo} . 


A simple example is the basin of attraction of a = 0 under the function 
f(z) = z 2 . This is the open unit disc because 

/"(z 0 ) —► 0 as n —► oo, for |z 0 | < 1, 
but 

f n (z o) -rtOasn-too, for |z 0 | > 1 (since |/"(z 0 )| > 1, 

for n = 1,2,...). 

Later in the unit we see some more complicated examples. 
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Some texts use 

(a) attractive or stable ; 

(b) repulsive or unstable; 

(c) neutral. 


Note that (d) is a special case of (a). 


Here, and subsequently, we may 
use z, rather than zo, as an 
initial term, when we do not 
need to label the sequence {z n }- 


Recall from Example 1.3(b) that 
/ n (z 0 ) = z\ , for n =0,1,2,... . 







Problem 1.7 


For each of the following functions /, determine the basin of attraction of the 
given fixed point a. 


(a) f(z) = ±z, a = 0 (b) f(z) = z 3 ,a = 0. 


1.3 Conjugate iteration sequences 

Consider the iteration sequence 

z n+ i = z 2 + 2z„, n = 0,1,2,..., with z 0 = — §, (1.3) 

and suppose that we wish to find a formula for z n in terms of n. This iteration 
sequence is quite complicated, and so it is sensible to look at the terms of {z n }: 

z 0 — — 5 > 2 1 ~ 4 i 2 2 — 16 > 2 3 — 256 ’ ' ‘ ‘ ' 

These terms suggest that z n —► — 1 as n —> oo, and so we make the substitution 
z n = w n - 1, 

and try to find a formula for w n . Substituting for z n and z n +i in (1.3), we 
obtain 

w n +i -l = (w n - l) 2 + 2 (w n - 1) 

= id 2 - 2 w n + 1 + 2 w n - 2 
= wl - 1 ; 

hence 

iu n+ i = w 2 . 

The iteration sequence 

w n +i — iy 2 , n = 0,1,2,..., with w 0 = z 0 + 1 = |, (1.4) 

is simpler than one given in (1.3) and, moreover, we know that 
Wn = (I) 2 > for n = 0,1,2,... . 

We deduce that the formula for z n in terms of n is 
z n =(|) 2 "- 1 , for n = 0,1,2,... . 

More generally, suppose that 

z n+1 = f(z n ), n = 0,1,2,..., 

is a given, but complicated, iteration sequence that we wish to investigate, and 
that h is a one-one function. Then 
w n =h(z n ), n = 0,1,2,..., 

is also an iteration sequence. Indeed, since z n = h~ 1 (w n ), we have 
w n+1 = h(z n+1 ) 

= h(f(z n )) 

= h(/(h _1 (u; n ))), for n = 0,1,2,... . 

Thus 

w n +\ = g(w n ), for n = 0,1,2,..., 
where the function g is given by g = h o / o h~ x (see Figure 1.6). With a 
suitable choice of h, the function g will be simpler than / (as in the above 
example, where f(z) = z 2 + 2z, g(w) = w 2 and h(z) = z + 1). 


Here z„ = / n (zo), where 
f(z) = z 2 + 2z. 


Thus w n = h(z n ), where 
w = h(z) = z + 1. 


Here io„+i = g(w „), where 
g(w) = w 2 . 

Example 1.3(b) 


The inverse function h 1 exists 
because h is one-one. 



h~ l 


Figure 1.6 g(w n )=(hofoh *)"(wn) 


We make the following definition. 




Definition The functions / and g are conjugate to each other if 
g = ho f o h -1 , 

for some one-one function h called the conjugating function. If the 
sequence {z n } is defined by 

*»+i =/(*«), n = 0,1,2,..., 

for some z 0 , and w n = h(z n ), for n = 0,1,2,..., then the sequence {w n } 
satisfies 

w n+1 = g(w n ), for n = 0, 1 , 2 ,..., 
and {z n } and {w n } are called conjugate iteration sequences. 


Remark In practice, the function g is usually found by substituting 

z n = h~ 1 (w n ) and z n+1 = h~ 1 (w n+1 ) into z n+1 = f(z n ) and then rearranging, 

as we did with (1.3). 

Since the sequence {m n } is the image of the sequence { z n } under the 
function h, it follows that both sequences have the same behaviour. For 
example, if both h and h~ l are continuous functions, then {z n } is convergent if 
and only if {w n } is convergent. 

Also, a is a fixed point of / if and only if h(a) is a fixed point of g. 

If the conjugating function h is required to be one-one and entire, then it must 
in fact be of the form h(z) = az + b, where a ^ 0. (This can be proved using 
the Casorati-Weierstrass Theorem.) When dealing with the iteration of 
rational functions (as in the next subsection), we also use extended Mobius 
transformations as conjugating functions, since these are one-one on C. 

Problem 1.8_ 

Prove that the iteration sequence 

z n+i=z n -zl, n = 0,1,2,..., 

is conjugate to the iteration sequence 

io„ +1 = n = 0,1,2,..., 

with conjugating function h(z) = -z -f If z 0 - |, what is w 0 ? 

Problem 1.9_ 

(a) Prove that the iteration sequence 

z n+ i=az n + b, n = 0,1,2,..., 
where a ^ 1, is conjugate to the iteration sequence 
w n+ i = aw n , n = 0,1,2,..., 
with conjugating function h(z) — z + b/(a - 1). 

(b) Hence obtain a formula for z n in terms of z 0 and describe the behaviour of 
the sequence {z n } when 

(i) |a| < 1; (ii) |a| = l,a / 1; (iii) |a| > 1. 


Do not confuse this use of the 
word ‘conjugate’, which is 
borrowed from group theory, 
with ‘complex conjugate’. 


Substituting w n = h(z n ), for 
n = 0,1,2,..., where h is a 
one-one function, amounts to 
making a change of variable. 


You dealt with the case a = 1 
Problem 1.3(a). 
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1.4 The Newton-Raphson method 


We shall now see how the ideas introduced so far help with the 
Newton-Raphson method, described in the Introduction. If p is a polynomial 
function, then the corresponding Newton-Raphson iteration sequence is 


n = 0,1,2,... 


, +1= , 

n+1 " A*nY 

We call the function to be iterated here the Newton-Raphson function 
corresponding to p, and denote it by N\ thus 


In general, N is a rational function (unless p is of degree 1) and has 
extension NtoC = Cu{oo}. Thus a Newton-Raphson iteration sequence may 
include the point at infinity among its terms. 


A key observation is that if a is a simple zero of p, then p'(a) ± 0 and 
N(a) = a. Thus a is a fixed point of N. To classify it, we evaluate N'(a): 


JV'(a) = 1 - 


jp'(a)) 2 -p(a)p"{a) 


(p'(a)) 2 


’ (p'(a)) 2 


(since p(a) = 0) 


Thus a simple zero a of p is a super-attracting fixed point for the 
Newton-Raphson function N. This is good news because it means that the 
Newton-Raphson method always converges rapidly to a simple zero a of p 
provided that our initial guess is close enough to a. 


In 1879, Cayley analysed the Newton-Raphson method when p is the quadratic 
function 

p(z) = z 2 + az + b, where a, b € C. 

If the function p has distinct zeros at a and f3, say, then these zeros must be 
simple. Thus a and (3 are super-attracting fixed points of the Newton-Raphson 
function N. By Theorem 1.1, there are open discs around a and (3 whose 
points are attracted to a and /3, respectively, under iteration of N. Cayley 
wished to know which points in C are attracted to a under iteration of N, and 
which are attracted to (3. In other words, what are the basins of attraction of a 
and (3 under N? 


Cayley found that the answer is remarkably simple: the perpendicular bisector 
of the line segment joining a and /3 forms a ‘watershed’ for this iteration 
process (see Figure 1.7). If zq falls on the same side of the watershed as a, then 
N n (z 0 ) —» a as n —» oo, but if z 0 falls on the other side, then N n (z 0 ) —* /3 as 
n — » oo. If zq falls exactly on the watershed line, then the sequence {N n (2 0 )} 
remains on the line! 

If we consider the Newton-Raphson function for p(z) = z 2 + az + b, 


N(z) = 


P(*) 

p>(z) 

z 2 + az + b 
2 z + a 


z 2 -b 
2 z + a’ 


then it is not at all evident why Cayley’s result should be true. However, we 
obtain a much simpler iteration sequence by using the conjugating function 



which is a Mobius transformation. The extended Mobius transformation h is 
one-one on C and maps q to 0 and /3 to oo. Also h has the remarkable property 


Unit Dl, Section 1 
Unit B3 , Section 5 


If a is a zero of order k, where 
k > 1, then it can be shown that 
a is a fixed point of N such that 
N'(a) = 1 — 1/fc < 1, 
so the Newton-Raphson method 
still works, but not so well. See, 
for example, Problem 1.11. 
Arthur Cayley (1821-1895) is 
well known for pioneering work 
in group theory and linear 
algebra. 



Figure 1.7 


Unit Dl, Section 2 
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that 


h(N(z)) = (, h(z )) 2 , 

which we ask you to prove in Problem 1.10. 

Thus, if 

z n+1 =N(z n ), n = 0,1,2,..., 
and if 

w n = h(z n ), for n = 0,1,2,..., 

then 

Wn+l = h(z n+1 ) 

= *(*(*»)) 

= ( h(z n )) 2 (by Equation (1.5)) 

= Wn, for n = 0,1,2,... . 

Thus 

w n = g n (w o), for n = 1,2,..., 

where g{w) = u; 2 . 

Now we know that 

—► 0 as n —> oo, if |w 0 | < 1, 

that 

w n —► oo as n —> oo, if |io 0 | > 1, 
and also that 

the sequence {w n } remains on the circle {w : |w| = 1} if |io 0 | = 1. 
To find out what happens to the original sequence {z n }, we note that 
h~ 1 (w n ) = z n , h~ 1 ( 0)=a, h~ 1 (oo)=P 

and 

ko! = l^o)| = |^^|. 

I 20 - p I 

Therefore we deduce from (1.6) and (1.7) that 
z„ -+ a as n -► oo, if |z 0 - a| < \zo ~ P\, 

that 

z n —» (5 as n —> oo, if \z 0 - a| > | z 0 - /?|, 
and also that 

{z n } remains on the extended line {z :\z — a\ = \z — /?|} U {oo} 
if zo is on this line. 

This is Cayley’s remarkable solution. 


(1.5) Identity (1.5) may be written 
the form 

N = h~ 1 ogoh, 
where g(w) = w 2 . 


(1.6) Example 1.3(b) 
and 

Problem 1.4(d) 

(1.7) 


Note that the function 
= ? 

is analytic at 0 and has a 
removable singularity at oo. 


Problem 1.10__ 

Prove Identity (1.5). 

{Hint: Note that both z = a and z = (3 satisfy the equation z 2 = -(az + b).) 

Problem 1.11_ 

Describe what happens under iteration of N if a = (3. 


Finally, we look briefly at the Newton-Raphson method for cubic polynomial 
functions. Here, you might guess that the complex plane divides itself into 
three simple regions, each surrounding one zero of p and consisting of those 
points which are attracted to that zero under iteration of the Newton-Raphson 
function N. Such indeed seems to have been Cayley’s hunch in 1879, although 
he was unable to prove such a result. With the help of computer-generated 
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pictures we can now see why Cayley had no chance of finding a simple solution 
to this problem! 

Consider p(z) = z 3 - 1, whose zeros are 

c*i =1, Q2 = e 2m ^ 3 = |(-1 + iV 3), a3 = e 4m ^ 3 = §(-l - iV 3). 

In this case 

w — ^ 

_ 2z 3 +1 
3z 2 

and, as before, aj, a 2 , Q3 are each super-attracting fixed points of N. 

Figure 1.8 shows in white the basin of attraction of the fixed point Qi = 1. 



Figure 1.8 The basin of attraction of 1 

The basins of attraction of ai, a 2 , 03 must be congruent to each other under 
rotation about 0 through 2n/3, because of the symmetry of <*1, a 2 , 03. But 
these basins are not at all simple, and they are not even regions (because they 
are, in fact, not connected). The union of these three strange basins is almost 
the whole of C. In addition, there is a complicated ‘watershed’ which separates 
the basins of attraction (see Figure 1.9) and which manages, somehow, to be 
the boundary of all three sets simultaneously! 



Figure 1.9 The watershed 


Thus the iteration of even fairly simple rational functions can lead to very 

complicated behaviour. In the next section, we find that complicated behaviour We return to the 

can occur even for the iteration of simple polynomial functions. Newton-Raphson method in 

Section 5. 
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2 ITERATING COMPLEX QUADRATICS 

After working through this section, you should be able to: 

(a) conjugate a given quadratic iteration sequence to one determined by a 
member of the family of functions {P c : c € C}, where P c (z) = z 2 +c; 

(b) explain why P c n (z 0 ) -► oo as n -► oo if \z 0 \ is large enough; 

(c) describe some of the properties of the escape set E c of P c and of its 
complement, the keep set K c ; 

(d) find the periodic points of P c , for certain values of c, and find their nature; 

(e) understand the definition of the Julia set J c . 

2.1 The basic quadratic family 

In Section 1 we saw that iteration sequences of the form 

z n+ i=az n + b, n = 0,1,2,..., (2.1) 

can be completely analysed; that is, for all a, 6 € C, we can describe the 
behaviour of {z n } for any initial term z 0 . In this section we begin to study 
iteration sequences of the form 

z„+i = az 2 + bz n + c, n = 0,1,2,..., (2.2) 

where a 0. We shall find that the family of such sequences is much more 
diverse than that given by (2.1). To begin with we note that every iteration 
sequence of the form (2.2) is in fact conjugate to one of a simpler type. 


Theorem 2.1 The iteration sequence 


z n+ i = az 2 + bz n + c, n = 0 , 1,2, ..., 


where a ^ 0, is conjugate to the iteration sequence 


Wn+i =wl + d, n = 0, 1,2, ..., 

(2.3) 

where d — ac + |6 — |6 2 . The conjugating function is 


h(z) = az+±b. 



Proof The recurrence relation (2.2) can be rearranged as 
az n+ i = ( az n + |6) 2 +ac- \b 2 , n = 0,1,2,... . 

Thus, putting w n = h(z n ), for n = 0,1,2,..., where /i(z) = az + ±6, we obtain 
w n+ i-\b = w 2 n +ac-\b 2 , n = 0,1,2,...; 
that is, 

tu n +i =w 2 +d, n = 0,1,2,..., 
where d = ac+\b - \b 2 , as required. ■ 

There are many different iteration sequences of the form (2.2) which are 
conjugate to any one iteration sequence of the form (2.3). This is illustrated in 
the following problem. 

Problem 2.1_ 

Use Theorem 2.1 to show that each of the following iteration sequences 

(a) z n+1 = 4z n (l - z n ), n = 0,1,2,..., with z 0 = §, 

(b) z n+1 = 1 - 2 z 2 , n as 0,1,2,..., with z 0 = 0, 
is conjugate to the iteration sequence 

w n+ i = w 2 - 2, n = 0,1,2,..., with wq = 0. 



Example 1.3(c) and Problem 1.9 


Problem 1.8 is a special case of 
this result, with a = -1, 6=1, 
c = 0. 


Multiply by a and complete the 
square. 
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Theorem 2.1 tells us that if we wish to understand the possible behaviour of 
quadratic functions under iteration, then it is sufficient to consider only those 
of the form (2.3), and it is convenient to relabel these as 

z n +i = z 2 +c, n = 0,1,2,..., 

where c is a complex parameter. We shall devote most of the rest of this unit 
to such iteration sequences, and so we introduce a name for the corresponding 
quadratic functions. 


Definition The set of functions {P c :cgC} defined by 
P c (z) = z 2 + c, 

where c € C, is the family of basic quadratic functions. 

For example, P 0 (z) = z 2 , P\(z) = z 2 + 1 and Pi(z) = z 2 + i axe all basic 
quadratic functions. In the following problems we ask you to establish some 
elementary properties of the basic quadratic functions. 

Problem 2.2_ 

(a) Determine the rules for the functions P 2 and P 2 . 

(b) Write down a formula for P" +1 (z) in terms of P”(z), and hence prove 
that P" is an even polynomial function of degree 2". 

Problem 2.3_ 

(a) Show that the fixed points of P c are \ ± yj\— c, and prove that at least 
one of these is repelling (unless c = |). 

(Hint: If a and (3 are the fixed points of P c , then 
±(P'(a) + P'(/3)) = l.) 

(b) What happens if c = |? 


2.2 The escape set and the keep set 

In Example 1.2(b) we observed that the iteration sequence 
z n = P£(zo), n = 0,1,2,..., 
can be determined explicitly as 

z n = zf, for n = 0,1,2,... . 

Thus 

z„ -* 0 as n -> oo, if |z 0 | < 1, ) 

z n -+ oo as n —► oo, if jz 0 | > 1, > (2.4) 

|z„| = 1, for n = 1,2,..., if |z 0 j = 1. J 

It is natural to ask whether a similar behaviour occurs for other values of c. We 
shall see later that when the initial term z 0 is small, the sequence 
z n = P?(z o), n = 1 , 2 ,..., behaves in dramatically different ways for different 
values of c. However, when the initial term z 0 is large these sequences behave 
in essentially the same way, for all values of c, as we now show. 


Remember that P c n 
nth iterate of P c . 


Po(z) = z 2 , 
so that c = 0. 





Theorem 2.2 Let r c = ± + yj\ + |c|. Then, for |z 0 | > nr, 
{|P”(2o)|} is an increasing sequence, 

and 

P?(zo) —> oo as n —> oo. 


Proof First note that, by the backwards form of the Triangle Inequality, 
|P c (z)| = | 2 2 + c|>|2| 2 -|c|. 

The number r c is the positive solution of the equation x 2 — |c| = x (see 
Figure 2.1), and we claim that if e > 0, then 

x 2 - \c\ > (1 + e)*, for x>r c + e. 

Indeed, if x > r c -f e, then 



x x 


>(r c + e)- 

_ r 2 + 2 r c e 


|c| 

r c + e 


W 


(2.5) 


(2.6) 


r c + 2 r c e + e 2 , . 2 i i \ 

=- r ' + e - (since r c — \c\ = r c ) 

r c +r c £ + e + e 2 t . N 

>- r + £ - (since r c > 1) 

= 1 + e, C 

as required for Inequahty (2.6). Inequalities (2.5) and (2.6) now give 
|Pe(*)|>(l+e)|*l>|*|, for \z\ >r c + £. 

If | 2 o| > r c + e, then we can apply this inequality successively to zq, P c (zo), 
P 2 (z 0 ),..., to deduce that the sequence {|P”( 2 0 )|} is increasing, and 

|P c n (2o)| > (l + e) n | 2 0 |, for n= 1,2,... . 


Hence P"(zo) —♦ oo as n —> oo. Since e is any positive number, the proof is 
complete. ■ 


Remarks 

1 Theorem 2.2 may be easier to understand if you think of oo as a fixed point 
of the extended function P c of P c . To discover the nature of this fixed point, we 
use the conjugating function h(z) = 1/2, which moves the fixed point from 00 
to 0. By this means the iteration sequence 
2 n+1 =2 2 +c, n = 0, 1 , 2 ,..., 


is conjugate to the iteration sequence {u>„}, where 

• (i-V— „=o,u,... 


As expected, because P c (00) = 00, the corresponding function 
Qc(w) = w 2 /(l + cw 2 ) has a fixed point at 0, which turns out to be 
super-attracting. Thus, by Theorem 1.1, 


w n = Qc( w o) —► 0 as n —» 00, provided that |w 0 | is small enough, 


and so 


2 n = P”(2 q) —> 00 as n —> 00, provided that |2 0 | is large enough. 


If c = 0, then r c = 1 and we 
recover the middle line of (2.4). 

y = x 2 -\c\ f 



Figure 2.1 r 2 - |c| = r c 


See Unit Dl, Subsection 2.2. 



Since Q' c {w) = - ———5, we 
(1 + cw 2 ) 

have <?c(0) = 0. 
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2 The formula r c = \ 4- yj\ + |c| will play a significant role later in the unit; 
therefore we indicate in Figure 2.2 the graph of the function |c| i—> r c , together 
with the graph of the identity function |c| i—> |c| for comparison. Notice, in 
particular, that 

|c| < r c ^ |c| < 2 (2.7) 


Problem 2.4_ 

(a) Verify that 

r 0 = 1, rj = |(l + \/5) and r_ 2 = 2. 

(b) Show that the sequences {P^(l)} and |P” 2 (2)} are constant. What does 
this tell you about the values r 0 and r_ 2 in relation to Theorem 2.2? 


We now investigate the set of all points which are attracted to oo or ‘escape’ 
to oo under iteration of P c . We call this set the escape set; its complement is 
the set of points which we ‘keep’ under iteration of P c . 

Definition For ceC, the escape set E c is 
E c = {z : P”(z) —> oo as n —> oo} . 

The keep set K c is the complement of E c . 

Remark The escape set E c may be thought of as the basin of attraction of oo. 
For example, we know by (2.4) that 

E 0 - jz : Po(z) = z 2 "-*ooasn-t oo j 
= {z:\z\> 1 }, 
and hence that 
K 0 = C - E 0 
= {z:\z\<l} 

(see Figure 2.3). 

This example is rather misleading as, for most values of c, E c does not have a 
simple shape. Figure 2.4 shows several examples of sets E c (in white!) and K c , 
plotted by computer. 



Figure 2.4 

In fact, c = —2 is the only other value for which E c and K c have a simple 
shape. We now ask you to investigate this case. 



Figure 2.3 

Here, and throughout this 
section, the square represented 
in such figures is 

{z : -2 < Rez.Imz < 2}. 




(c) Ec and K c :c=-1 
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L is the interval [—2,2], 


Problem 2.5_ 

Let L be the line segment {x + iy: |x| <2 ,y = 0}. 

(a) Let 

z n+1 =zl- 2 , n = 0,1,2,.... 

Prove that 

if z 0 € L, then z n € L, for n = 1,2,...; 

if z 0 € C - L, then z n 6 C — L, for n = 1,2,... . 

(b) Prove that if z 0 G C - L, then the sequence {z n } in part (a) is conjugate to 
an iteration sequence of the form 

w n+ i = wl, n = 0,1,2,..., 

and deduce that z n —» oo as n —> oo. 

(Hint: Recall from Unit Dl, Subsection 4.4, that the Joukowski function 
J(w) = w + l/w maps {to : |to| > 1} one-one and conformally onto C - L, 
put w n = J~ l (z n ), for n = 0,1,2,..., and verify that J(to n+1 ) = J(to£).) 

(c) Deduce that £_ 2 = C - L and R_ 2 = L (see Figure 2.5). 


Though the set E c is usually very complicated, a number of general 
observations can be made about it. For example, we can show that E c has the 
property of being completely invariant under P c . 


Definition A set A is completely invariant under a function / if 

z e A f(z) e A. 


This means that if z lies in A, then f(z) and all its iterates also lie in A and, 
moreover, any point whose iterates eventually lie in A must itself he in A. For 
example, it is easy to check that each of the sets C, {0}, C - {0}, {z : \z\ < 1}, 
{z : |z| = 1} and {z : |z| > 1} is completely invariant under the function 
/(*) = * 2 - 

The following result lists several key facts about E c and K c . 


Theorem 2.3 For each c € C, the escape set E c and the keep set K c have 
the following properties: 

(a) E c D{z: |z| > r c } and K c C {z : |z| < r c }; 

(b) E c is open and K c is closed; 

(c) E c jz C and K c 0; 

(d) E c and K c are each completely invariant under P c - 

(e) E c and K c are each symmetric under rotation by n about 0; 

(f) E c is connected and K c has no holes in it. 


Remarks 

1 In each case the property of K c is equivalent to the corresponding property 
of E c . Thus in the proof we need only establish the results for E c . 

2 Note that Properties (a), (b) and (c) combine to tell us that K c is a 
non-empty compact set. 


E. 2 

_ A'-; [ 

-2 2 


Figure 2.5 


The statement that K c has no 
holes in it just means that E c is 
connected. 
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Proof of Theorem 2.3 

(a) We have already proved part (a) in Theorem 2.2. 

(b) Suppose that z 0 S E c . Then, by definition, P"(z 0 ) —» oo as n —► oo. Thus, 
for some n 0 , |P"°(zo)| > r c . 

Let e = |P c n °(z 0 )| - r c . Since e > 0 and P”° is a polynomial function, P c n ° 
is continuous at zo and so there exists 6 > 0 such that 

\z-z 0 \<6 =* \P?°{z) - P c n °(z 0 )| < £ 

and hence 

|*-*>|<« =* |P”°(z)|>r c 

(see Figure 2.6). 

It follows from Theorem 2.2 that 

|z — z 0 1 < 5 => P"(z) —* oo as n —> oo, 

so that {z : |z — zo| < 5} C E c . This shows that E c is open. 

(c) The set E c is not the whole of C, because it does not contain the fixed 
points of P c , found in Problem 2.3. (These must lie in K c .) 

(d) We ask you to prove part (d) in Problem 2.6. 

(e) Since 

P?(-z) = P c n (z), for n = 1,2,..., 
because P™ is an even function (see Problem 2.2(b)), we have 
P c "(-z) —► oo <==> P”(z) —+ oo 
and so 

-z € E c z e E c . 

(f) Finally, we show that E c is connected, that is, each pair of points in E c can 
be joined by a path lying in E c . Since the set A c = {z : \z\ > r c } is 
connected and A c C E c , it is sufficient to show that each point a in E c can 
be joined to some point of A c by a path in E c . We prove this by 
contradiction. 

Suppose in fact that a is a point of E c which cannot be joined to the 
set A c in this way. We define the following set: 

TZ = {z e E c : z can be joined to a by a path in E c } 

(see Figure 2.7). 

Then TZ ^ 0 (since a € TZ), TZ is open (because if z can be joined to a 
in E c , then so can points of any open disc in E c with centre z) and TZ is 
connected (because pairs of points in TZ can be joined in TZ via a). Thus TZ 
is a subregion of E c . Since a cannot be joined in TZ to A c and TZ is open, 
we deduce that TZ C {z : \z\ < r c }. 

We can now use the Maximum Principle. If /3 e dTZ, then 0 does not lie in 
E c (since, otherwise, we could enlarge TZ slightly). Thus 

\P?(P)\ < r c , for n= 1,2,... . 

By applying the Maximum Principle to each of the polynomial functions 
P c n on TZ, we obtain 

^ r c , for n — 1,2,..., and z € TZ, 
which contradicts the fact that TZ C E c . Hence E c is connected. ■ 

Problem 2.6___ 

Prove part (d) of Theorem 2.3. 


This proof may be omitted 
first reading. 


P"°(z) 



Figure 2.7 
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2.3 Periodic points 

Figure 2.4 shows that the keep sets K c have a remarkably diverse form. In 
order to investigate the shape of K c , we need to identify as many points in K c 
as possible. We already know that the fixed points ~ ± \J\ — c of P c he in K c , 
and we can find other points of K c by generalizing the notion of a fixed point. 
As an example, consider the function P-i{z) — z 2 - 1, whose two fixed points 
|(1 ± \/5) must lie in K-\. In addition, the points 0 and — i, have the 
property that 

P_i(0) = -1 and P_i(-1) = 0 (2.8) 

(see Figure 2.8). Thus the sequence P™ 1 (0),n = 0,1,2,..., cycles endlessly 
between the points 0 and -1, as does the sequence P”j(—1), n = 0,1,2,... . 
This means that 0 and -1 must both lie in K-\. Since K-i is symmetric 
under a rotation by 7r about 0, we can begin to build up a picture of K_ x (see 
Figure 2.9). This is very far from the complicated set in Figure 2.4(c), but at 
least it is a start! 


Problem 2.3(a) 



Figure 2.8 


Theorem 2.3(e) 


§(1 - \/ 5 ) ' 
, , 1 


-1 0 

|(1 -t- >/ 5 ) 


Figure 2.9 


The fact that 0 and —1 satisfy Equations (2.8) can be interpreted as saying 
that 0 and —1 are both fixed points of the second iterate 

P* 1 (z) = P. 1 (z 2 - 1) 

= (* 2 -!) 2 -l 

= z 4 - 2z 2 . 


Indeed, it is evident that P 2 j(0) = 0 and Pl 1 (-1) = -1. 

Points of this type, which are fixed points of higher iterates of a function /, are 
called periodic points of /; they repeat periodically under iteration of /. 


Definition The point a is a periodic point, with period p, of a 
function / if 

f p (oc) = a, but f k (a) ^ a, for k = 1 , 2 , ... ,p — 1. 

The p points 

a. /(<*)> / 2 (“). •■•1 / p-1 (<*) 

then form a cycle of period p, or a p-cycle of / (see Figure 2.10). 


Remarks 

1 Note that if we apply / repeatedly to points of a p-cycle, then we just obtain 
points of the p-cycle. Also the points of a p-cycle must be distinct. Indeed, if 
we had 


Thus p is the smallest positive 
integer such that 
/» = a. 


;/» 

•/V) 


f k (a) = f e (a), where 0<fc<£<p-l, 


Figure 2.10 


then the point f p (a) = f p k {f k {ot)) would have to lie among the terms 
f k+1 (a),..., f e (a), 


which do not include a, which is a contradiction. Therefore, the points of a 
p-cycle are each distinct periodic points of / with period p. 

2 Note that a fixed point of / is a 1-cycle of /. 

3 Evidently, all the periodic points of P c lie in the keep set K c . 
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Determining the periodic points with period p of a given function /, for a given 
p > 1, is usually more difficult than determining the fixed points of /. This is 
because the equation 

f p (z) = z (2.9) 

usually has many more solutions as p increases and the rule for f p is usually 
more complicated. Notice, however, that not all solutions of Equation (2.9) 
need be periodic points with period p. For example, any fixed point of / also 
satisfies Equation (2.9). More generally, if q is a factor of p, then any solution 
of f g {z) = z is also a solution of f p (z) = z. 


Example 2.1 

Determine all periodic points with period 2 of the function Po(z) = z 2 , and 
write down the corresponding 2-cycles. 

Solution 

Since Pq(z) = (z 2 ) 2 = z 4 , we have to solve the equation z 4 = z: 


*=*z(z 3 -l)=0. 

The solutions of this quartic equation are 0,1, e 2 "/ 3 = |(—1 + \/3i) and 
e 4 "/ 3 = |(—1 — y/3i). Of these, the points 0, 1 are fixed points of Po, whereas 

F’o^e 2 ’"/ 3 ) = (e 2 ”/ 3 ) 2 = e 4iri / 3 , 

and 

Po(e 4 " /3 ) = (e 4ir< / 3 ) 2 = e 2 -/ 3 . 

Hence both e 2 "/ 3 and e 4m / 3 are periodic points of P 0 with period 2, and they 
belong to the 2-cycle e 2 ”/ 3 , e 4nl / 3 . ■ 

Note that, as expected, the 2-cycle found in this example lies in the keep set 
Kq = {z: \z\ < 1} (see Figure 2.11). 



Figure 2.11 


Problem 2.7_ 

(a) Prove that —i is a periodic point with period 2 of the function 

Pi(z) = z 2 + i. Hence find five points in Ki, none of them fixed points of The set Ki was plotted in 
Pi, and plot them. Figure 2.4(b). 

(b) Determine all periodic points with period 3 of the function P 0 (z) = z 2 , and 
write down the corresponding 3-cycles. Plot all the fixed points, 2-cycles 
and 3-cycles of P 0 on the same diagram. 

(c) Prove that |(—1 + y/2) is a periodic point of P_ 5 / 4 . 


We now return to the function P-i(z) = z 2 - 1, which has the 2-cycle 0,-1; 
that is, 0 and —1 are fixed points of the second iterate 

Pi 1 {z) = z 4 -2z 2 . 

Since 

(P! 1 )'(z)=4z 3 -4z 

= 4z(z 2 — l), 

it follows that 

(P 2 1 )'(0) = 0 and (P! 1 )'(-1) = 0, 

so both 0 and -1 are super-attracting fixed points of P 2 X . The fact that these 
two derivatives have the same value is no accident^ as the following result 
shows. 
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Theorem 2.4 Let a, f(a), / 2 (a), ■. . ,/ p 1 (a) form a p-cycle of an 
analytic function /. Then 

(a) ( f p )' (a) = f'(a) x /'(/(a)) x f'(f 2 (a)) x • ■ • x /'(/^(a)) (2.10) 
and, moreover, 

(b) the derivative of f p takes the same value at each point of the p-cycle; 
that is, 

(/ p )' (<*) = (/ p )' (/(a)) = (/ p )'(/ 2 («)) = ... = (/ P )'(/ P_1 (<*)). 


Proof Since f p (z) = /(/(• • ■ (f(z)) ■ ■ •)), where the function / is applied p 

times, we deduce from repeated applications of the Chain Rule that Unit A4, Theorem 3.1 

{f P )' (*) = f'{F-\z)) x • ■ ■ x /'(/(*)) x f(z). 

By putting z = a, we obtain (a). 

Thus ( f p )' (a) is the product of the derivatives of / at the points of the p-cycle, 
and so ( f p )' (/(a)) is also the product of the derivatives of / at the points of 
the p-cycle, and similarly for the other points f 2 (a), ..., / p-1 (a). This 
establishes (b). ■ 

Theorem 2.4 allows us to classify the periodic points of an analytic function / 
by using the number ( f p )' (a), which is called the multiplier of the 
corresponding cycle. We shall see shortly that different types of cycle lie in 
different parts of the keep set. 


Some texts use the name 
eigenvalue rather than 
multiplier. 


Definitions If a is a periodic point, with period p, of an analytic 
function /, then a and the corresponding p-cycle are 

(a) attracting, if |(/ p ) , (a)| < 1; 

(b) repelling, if |(/ p )'(a)| > 1; 

(c) indifferent, if | (/ p )'(a) | = 1; 

(d) super-attracting, if ( f p )' (a) = 0. 

In the next example we demonstrate two ways to calculate the multiplier. 

Example 2.2 

Determine the nature of the periodic point e 2,ri / 3 of Po(z) = z 2 . 

Solution 

In Example 2.1 we found that e 2 ’ 1 ’*/ 3 is a periodic point, with period 2, of the 
function Pq. Since P§(z) = z 4 , the multiplier is 

(P 2 ) , (e 2 "/ 3 ) = 4(e 2 "/ 3 ) 3 (since (P 2 )' (z) = 4z 3 ) 

= 4, 

so that |(P 0 2 ) (e 2?rt / 3 )| > 1. Thus e 2ir ’/ 3 is a repelling periodic point of P 0 . 
Alternatively , the point e 2ni is part of the 2-cycle e 2ir */ 3 , e 4 "/ 3 for the 
function P 0 . Thus, by Theorem 2.4(a), the multiplier of the 2-cycle 

e 2iri/3 e 4iri/3 j g 

^o ( e2wi/3 ) Po(e 4iri/3 ) = fee 2 ™' 3 ') (2e 47ri / 3 ) (since P^z) = 2z), 

= 4, 

as before. ■ 

Note that the second method avoids the calculation of the rule for the pth 
iterate, which can be involved. 
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Problem 2.8- 

Determine the nature of each of the following periodic points, found in 
Problem 2.7. 

(a) —i, a periodic point of P* 

(b) e 2 "/ 7 , a periodic point of P 0 

(c) |(-1 + \/2), a periodic point of P- 5/4 


We now look more closely at where the periodic points of P c lie in K c . As an 
example, recall that the function P-i(z) — z 2 — 1 has fixed points at 
|(1 ± y/5), which are both repelling because 

|P: 1 (|(1 + V5))| =1 + V5>1 and |(P_ 1 )'(K 1 - >/5))| = V5 - 1 > 1 . PL 1 (z) = 2z 

Also, we saw before the proof of Theorem 2.4 that P_i has the super-attracting 
2-cycle 0, —1. In Figure 2.12, these points are plotted on a computer-generated 
picture of dK- X . 



This picture suggests that the repelling fixed points |(1 ± \/5) lie on the 
boundary of K- X , whereas the super-attracting 2-cycle 0, —1 lies in the interior 
of K-i- A similar phenomenon occurred for P 0 in the solution to 
Problem 2.7(b). The super-attracting fixed point 0 of P 0 lies in the interior 
of K 0 , whereas the repelling fixed point at 1 and the 2-cycle and the 3-cycles, 

which are all repelling, lie on the boundary of K 0 . In fact, the following general See Problem 2.8(b), for example, 
result holds. 


Theorem 2.5 Let a be a periodic point of the function P c . 

(a) If a is attracting, then a is an interior point of K c . 

(b) If a is repelling, then a is a boundary point of K c . 


The definitions of interior point 
and boundary point are given in 
Unit AS, Subsection 5.1. 


This proof may be omitted on a 

(a) Suppose that a is an attracting periodic point of P c , with period p. Then a ® rst readin 8- 
is an attracting fixed point of the pth iterate PJ\ Hence, by Theorem 1.1, 

there is an open disc with centre a whose points are attracted to a under 
iteration of Pf. These points therefore do not escape to 00 under iteration 
of P c , and so they must lie in K c . Hence a is an interior point of K c . 

(b) First suppose that a is a repelling fixed point, so that 

P c {a) = a and |P^a)| >1. 
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Since a € K c , we need to show that a is not an interior point of K c . If it Since K c is closed, 
were an interior point, then we could choose an open disc {z : \z — a\ < r} 8 K C = K c — int K c - 
lying in K c . In that case 

P"(z) e K c , for \z- a\= \r and n= 1,2,, 
by Theorem 2.3(d), and hence 

\P?( Z )\ < r c , for \z - a\ = \r and n = 1,2,..., r c = I + + |c| 

by Theorem 2.3(a). Now we apply Cauchy’s Estimate to each of the Unit B2, Problem 3.3 

polynomial functions P” to deduce that 

|(P”)'(a)| <y- = — , fern = 1,2,.... (2.11) 

2 r r 

By the Chain Rule, as in the proof of Theorem 2.4, 

(P?y (a) = P'(Pr») X • • ■ x P c '(P c (a)) x P'(a) 

= (*>»)”, 

since, by assumption, a is a fixed point of P c . Because |P'(a)| > 1, this 
implies that the sequence 

|(P")'(a)| = |P'(a)|", n = 1,2,..., 
tends to oo, contrary to Estimate (2.11). Thus a e dK c . 


If a is a repelling periodic point of P c , with period p, then a similar 
argument applies (with Pf rather than P c ); we omit the details. ■ 


Remarks 

1 Notice in part (a) that each point of the cycle 
a, P c (a), P c 2 (a),..., PrV), 


is an attracting periodic point of P c , with period p, by Theorem 2.4(b), and so 
each point of this cycle is an interior point of K c . The effect of P c is to map 
any point z near a to a point P c (z) near P c (a), and so on, round and round the 
cycle (see Figure 2.13, in which, for convenience, each disc has the same radius). 



In this way the sequence {P"(z)} forms itself into p convergent subsequences, 
each converging to a point of the attracting p-cycle. In fact it can be shown 
that every interior point of K c will be attracted in this way to the attracting 
p-cycle. See Figure 2.14, which shows an interior point zq of K-\ being 
attracted to the super-attracting 2-cycle 0,-1. 

2 In view of Theorem 2.5, it is natural to ask where in K c do any indifferent 
periodic points of P c lie. The answer is that it depends in a rather complicated 
way on the multiplier of the periodic point. This multiplier has modulus 1, and 
so it is of the form e 2ma , where 0 < a < 1. It can be shown that if a is a 


If the sequence 

z,P!(z),P?(z),... 
tends to a, then the sequence 
P c {z),P’ + 1 (z),P? p+ 1 {z),... 
tends to P c {oc), by the 
continuity of P c at a. 



Figure 2.14 

dK _! and z n = P2 1 {z 0 ) 
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rational number, then the periodic point lies on dK c , whereas if a is an 
irrational number, then the periodic point usually lies in the interior of K c . 

3 If we know that the set K c has no interior points, for a given value of c, then 

Theorem 2.5(a) tells us that P c can have no attracting periodic points. For 

example, K -2 = {x + iy : |x| < 2, y = 0} has no interior points and so Problem 2.5 

P_ 2 (z) — z 2 — 2 has no attracting periodic points, a fact which is not 

immediately obvious! 


2.4 The Julia set of P c 


We saw in the previous subsection that the boundary of the keep set K c (which 
is also the boundary of the escape set E c ) contains all the repelling periodic 
points of P c . This set is of particular interest because it forms the ‘watershed’ 
between those points which escape to 00 under iteration of P c and those which 
are kept (see Figure 2.15). It is called the Julia set of P c , in honour of the 
mathematician who first studied it in detail. 



(a) K c :c= -0.5 + 0.5i 



This section is intended for 
reading only. 


Gaston Julia (1893-1978), a 
French mathematician, 
developed the basic theory of 
the iteration of analytic 
functions while recovering from 
wounds received in World 
War I, and won a prize for this 
work from the Parisian 
Academy of Sciences. He was a 
professor at the Sorbonne, and 
worked on many aspects of 
complex analysis. 


Figure 2.15 


Definition The Julia set J c of P c is the boundary of K c . 


As noted earlier, the set K c has no holes in it and so K c can be thought of 
as J c together with the ‘inside’ of J c . Hence K c is often called the filled-in 
Julia set For example, 

K 0 = {z : \z\ < 1}, so J 0 = dK 0 = {z : \z\ = 1}, 

and 

K -2 = {x + iy : \x\ < 2, j/ = 0}, so J _2 = OK -2 = K- 2 - Notice that if K c has no interior 

A number of general properties of J c can be deduced from Theorem 2.3. For P oints then K c . 
each c e C, the Julia set J c is a non-empty compact subset of {z : |z| < r e }, 
which is completely invariant under P c and symmetric under rotation by n 
about 0. Of these properties, only the complete invariance is a little tricky to 
prove, and this holds because of the following facts: 

1 K c is completely invariant (Theorem 2.3(d)); 

2 int K c is completely invariant (by the continuity of P c and the Open Unit C2, Theorem 3.1 

Mapping Theorem); 

3 J c = dK c = K c - int K c . 
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The definition of J c could be used to plot J c , but there axe also other methods. 

We have already seen that J c contains all the repelling periodic points of P c , Theorem 2 . 5 (b) 

and so a knowledge of a number of these points gives us some information 
about the shape of J c . In fact it can be shown that 


J c is the smallest closed set which contains all the repelling 

periodic points of P c . (2.12) 

Thus the shape of J c is entirely determined by these points. However, 
calculating a large number of the repelling periodic points would be a difficult 
task in practice. 



A more satisfactory method of plotting J c is to use the complete invariance of 
this set under P c . This complete invariance tells us that if c*i € J c , then the 
solutions of the equation P c (z) = c*i, that is, z 2 + c = a 1: also lie in J c . This 
equation has the solutions ±y/a\ — c, which are two new points, c* 2 , 03 say, of 
J c (see Figure 2.16). Now, however, we can repeat this process with 0 : 2 , a 3 
instead of ot\ to obtain four new points 04 , a 5 , a 6 , a 7 in J c . This process, 
which is illustrated schematically in Figure 2.17, is known as backward 
iteration. It can be shown that 


Figure 2.16 

Exceptional cases occur when 
ai = c, or when ai is a fixed 
point of P c . In these cases, there 
is only one new point. 


J c is the smallest closed set which contains all the backward iterates 
of any given point of J c . (2-13) 

Thus the shape of J c is entirely determined by these backward iterates. The 
calculation of such backward iterates is not quite straightforward, even with a 
computer, because the tree-like structure shown in Figure 2.17 needs careful 
handling. A common short cut is to make a random choice of square root at 
each level and plot the resulting sequence; for example, 


..A... ...A. 

A A A A 


<*1, <*3, <*6, <*12, .... 

A convenient choice of starting point c*i for the backward iteration is a 
repelling or indifferent fixed point of P c . You can see the result of this method, 
for several values of c, in Figure 2.18. 



Figure 2.18 



(b) J c : c = -0.123 + 0.745i (rabbit) 



Remarks 

1 Property (2.12) implies that near every point of J c , we can find repe lling 
periodic points. The effect of these repelling periodic points is to make the 
behaviour of the function P c on J c extremely unstable, in the sense that points 
near together on J c tend to be pushed apart under the iteration of P c . Such 
behaviour is often described as chaotic. 

By contrast the behaviour of P c on C - J c is stable ; that is, points close 
together in C - J c behave in essentially the same way under iteration of P c . 
This distinction between stable and unstable behaviour can be used to define 
the- notion of a Julia set for any entire function or rational function. 
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2 Julia sets display a remarkable ‘self-similarity’ property. For each c, the 
shape of any part of J c appears to be repeated all over J c and is seen even 
when we zoom in closer and closer to J c . This is a consequence of the complete 
invariance of J c together with Property (2.13) of J c . As a result, Julia sets are 
often described as fractals. The name fractal was introduced by B. Mandelbrot 
(of whom, more later) in 1975 to describe a type of set which is extremely 
irregular, and yet has an underlying structure that can be seen under 
successive magnifications of the set. The exact definition of a ‘fractal’ is still 
under discussion, but it has to do with certain methods of measuring the 
dimension of a set which may give non-integer answers! 


3 GRAPHICAL ITERATION 


After working through this section, you should be able to: 

(a) use graphical iteration to determine the behaviour of real iteration 
sequences; 

(b) describe properties of the keep sets K c , for ceR, which can be obtained 
by using graphical iteration. 

In this section we make some observations about the nature of the keep sets K c 
when c is a real number. One simple observation is that if c is real, then K c is 
symmetric under reflection in the real axis (see, for example, K\ in Figure 2.4). 
This holds because, for ceR, 

P c (z) = z 2 +c=z*T2=P&), for zee, 
so that, for c S IR, 

P”(z) —> oo as n —> oo Pc{z) —» oo as n —> oo; 

thus, z € K c if and only if z € K c . 

We obtain some more interesting results by using a technique called graphical 
iteration, which applies only to the iteration of real functions. 

3.1 What is graphical iteration? 

If / is a real function, then an iteration sequence of the form x n+1 = f(x n ) can 
be represented graphically by using the two graphs y = f(x) and y — x plotted 
together, as in Figure 3.1. Note that any point where y = f(x) meets y = x 
corresponds to a fixed point of / (for example, the point a in Figure 3.1). 



These two diagrams illustrate a two-stage process for finding geometrically the 
position on the x-axis of the term x n+ i, given the position of the term x n : 

(a) draw a vertical fine to meet y = f(x) at (x n ,/(x n )) = (x n ,x n+1 ); 

(b) draw a horizontal line to meet y = x at (x n+ i,x n +i). 






Given any initial term Xo, we can apply the above process repeatedly to 
construct the sequence {x n } geometrically, and thus obtain information about 
the behaviour of {x„}. For example, with the function / in Figure 3.1 and with 
xo = 0, we obtain the behaviour illustrated in Figure 3.2, which strongly 
suggests that {x n } tends to the fixed point a. 

In the following example we carry out graphical iteration with a particular 
function /. 

Example 3.1 

Let /(x) = |x + 1. 

(a) Plot y = f(x) and y = x on the same diagram and use graphical iteration 
to plot the iteration sequences 

x n +i — /(x n ), n = 0,1,2,..., 
with xo = 0 and xo = 3. 

(b) Describe the behaviour of each of these sequences {x n } and check that 
your answer agrees with the solution to Problem 1.9(b)(i). 

Solution 

(a) The graphs y — f(x) and y = x are plotted in Figure 3.3. They meet at the 
point (2,2), which corresponds to 2, the only fixed point of /. The 
sequences x n+1 = /(x n ), n = 0,1,2,..., with x 0 = 0 and x 0 = 3, are also 
plotted. 



(b) For both values of xo, the figure strongly suggests that x n —> 2 as n —> oo. 
In Problem 1.9(b)(i) we found that if |a| < 1, then the iteration sequence 
z n+ \ = az n + 6, n = 0,1,2,..., converges to (the fixed point) a = 6/(1 — a) 
for all z 0 . Here we have a = |, b = 1 and a = 2, so our answer agrees with 
this result. ■ 

The next problem gives you a chance to try out graphical iteration. 

Problem 3.1_ 

Let /(x) = — 2x + 1. 

(a) Plot y = /(x) and y = x on the same diagram and use graphical iteration 
to plot the iteration sequences 

x„+i = f(x„), n = 0,1,2,..., 

with x 0 = 0 and x 0 = 3. 

(b) Describe the behaviour of each of these sequences {x n } and check that 
your answer agrees with the solution to Problem 1.9(b)(i). 



Figure 3.2 
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3.2 Real quadratic iteration sequences 

We now apply graphical iteration to real quadratic iteration sequences of the 
form 

X n+1 = P c (x n ) = X 2 n + C, 71 = 0,1,2,..., 

where c is real. A great deal can be said about sequences of this special form, 
but we confine attention to a few basic results which throw some light on the 
corresponding keep sets K c . 

To illustrate the method, we consider the case c = 0. The graphs 
y = P 0 (x) = x 2 and y = x meet at (0,0) and (1,1), corresponding to the fixed 
* points 0, 1 of P 0 , and the iterations plotted in Figure 3.4 indicate that if 

|x| < 1 then 

0 < Po(x) <1, for n = 1,2,... . 

On the other hand, if |x| > 1, then 
P” (x) - oo as n - oo. 

These results show that the part of the keep set K 0 = {z : \z\ < 1} on the real 
axis is the closed interval [-1,1], as expected. 

Now the only values of c for which we know K c explicitly are c = 0 and c = -2. 
For other real values of c, we may expect that graphical iteration will give new 
information about K c . To see what kind of information can be obtained in this 
way, try the following problem. 

Problem 3.2- 

(a) Plot y = x 2 + 1 and y = x on the same diagram and apply graphical 

iteration to the sequence x n+ i = + 1 , n = 0, 1 , 2 ,..., with your own 

choice of initial term x 0 . 

(b) Explain why the sequence {x n } tends to infinity. 

(c) What do you deduce about the set K\1 


The solution to Problem 3.2 suggests that the presence or absence of real fixed 
points of P c {x) = x 2 + c is of fundamental importance to the behaviour of 
iteration sequences of the form x n+1 = x\ + c. Since the fixed points of P c are 
1 ± \J\ ~ c, the following lemma is evident. 


Lemma 3.1 If c € R, then the function P c has 

(a) no real fixed points if c > 

(b) the single fixed point |, if c = 

(c) the two real fixed points | ± \J\ — c, if c < 


Part (a) of Lemma 3.1 shows that if c > then the graph y = x 2 4- c lies 
entirely above y = x (see Figure 3.5). It follows by graphical iteration that if 
c > |, then 

P"(x) -* oo as n -mx, for all x € R. (3.1) 

Thus, if c > then all real values of x escape to oo under iteration of P c , and 
so no real values of x belong to the keep set K c . This gives the following result. 


Theorem 3.1 If c > then K c n R = 0. 



Figure 3.4 


Problem 2.3 



Figure 3.5 c > j 


An algebraic proof of (3.1) is as 
follows. 

If x n +i = Xn + \ + £, where 
e > 0, then 

x n +i -x n = (x„ - \) 2 + e > e; 
hence 

x n > Xo + ne 
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It follows, by Theorem 2.3, parts (c) and (e), that if c > ±, then K c must have 
points in both the upper and lower open half-planes. Thus, for these values of 
c, the set K c is in at least two pieces, and so is not connected. Moreover, the 
plot of Ki in Figure 2.4 suggests that the set K l does not meet the imaginary 
axis. We ask you to prove this for K c ,c > ±, in the next problem. 

Problem 3.3___ 

(a) Show that if c is real and y is real, then P c (iy) is real. 

(b) Deduce from part (a) and Theorem 3.1 that if c > ±, then K c does not 
meet the imaginary axis. 


If c < 11 then, by Lemma 3.1, P c has either one or two real fixed points and so 
the keep set K c does meet the real axis. It turns out that if c lies in the interval 
[ — 2, |], then the set K c n R is precisely equal to the symmetric closed interval 



Theorem 3.2 if -2 < c < i, then K c n R = i c . 


Proof The graphs y = P c {x) and y = x are shown in Figure 3.6, together with 
a square S with sides parallel to the axes, which meets both axes in the 
interval I c . The key to the proof is the observation that if -2 < c < then the 
points of the graph y = P c (x), for x € I c , lie in S. If 0 < c < |, then this is 
evident because y — P c (x) does not extend below the x-axis. If -2 < c < 0, 
then we need to show that the lowest point of y = P c (x) does not lie below the 
bottom edge of the square 5, that is, 

c> - \J\~c, for —2 < c < 0. (3.2) 

This inequality can be verified directly or we can use Figure 2.2. This shows 
that 


Note that I 0 = [-1,1] and 
1-2 = [—2,2], as expected. 




and hence that 


for 0 < |c| < 2, 


which gives Inequality (3.2). 


It follows that 


x S I c => P c n (x) e 7 C , for n = 1,2,..., 


! 

I 


Moreover, graphical iteration shows that 
x £ I c => Pc(x) ->ooasn-*oo 
=*x(£K c . 

Thus K c nU = I c . U 

If c = -2, then we know that 

K - 2 = {x + iy : |x| < 2,y = 0} = /_ 2 ._ 

However, it can be shown that for other values of c in the interval [-2, |], the 
keep set K c does not lie entirely on the real axis. 
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For c < -2, the keep set K c has a more interesting intersection with the real 
axis, which is described in the next result. 


Theorem 3.3 If c < —2, then the set JT c n[R consists of the closed 
interval I c from which a sequence of disjoint, non-empty, open 
subintervals of I c has been removed. In particular, 0 £ K c . 


Proof Let S be the square used in the proof of Theorem 3.2, which meets the 
axes in the interval I c . First note that, by graphical iteration, 

P?(x) -* oo as n —> oo, for x € IR - I c , (3.3) 

so that K c n IR C I c . 

Now, since c < —2, the lowest point ony = P c (x) lies below S, as in Figure 3.7. 
Therefore, the set Ao of points in I c which escape from I c after exactly one 
iteration of P c , 

A 0 = {xCl c : P c (x) i I c }, 

is an open subinterval of I c with centre 0 (see Figure 3.7). In view of (3.3), it 
follows that 

P"(x) -» oo as n -too, for a: € .Ao, 
and so the points of Aq do not lie in K c . In particular, 0 £ K c . 



Figure 3.7 



Figure 3.8 


We now consider the set A\ of points in I c which remain in I c for exactly one 
iteration of P c : 

A 1 = {x€l c : P c (x) G I c , but P c 2 (x) i I c ) 

= {xel c : P e (x) G A 0 }. 

The set A\ (see Figure 3.8) consists of two open subintervals of I c which are 
positioned symmetrically on either side of 0 and contain the points ±y/—c, 
which are the zeros of P c . 

More generally, we consider the set A n of points in I c which remain in I c for 
exactly n iterations of P c , defined inductively as follows: 

A n = {x G Ic : P e {x) € A„_i}, for n = 1,T,... . 


The construction lines with 
single and double arrowheads in 
Figure 3.8 indicate how the 
endpoints of the two open 
subintervals, which comprise Ai, 
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Since each open interval in A n _i gives rise to two open intervals in A n , it This can be proved by using 

follows that the set A n consists of 2" disjoint open subintervals of I c . Now, any Mathematical Induction, 
point of I c which escapes to oo under iteration of P c must lie in exactly one of 
the sets A n . Thus the sets A n are disjoint and 

K c D IR = I c - (A 0 U A! U • • •), 
which gives the required structure. ■ 

Remarks 

1 The set K c n R = I c - (A 0 U Ai U • • ■) is infinite because it contains all the 
endpoints of all the intervals which comprise A n , for n = 0,1,2... . Actually, 
there are infinitely many points in K c D R which are not endpoints of this type, 
but these are harder to identify. 

2 It can be shown that if c < —2, then K c C IR (as is the case for c = —2). 

Problem 3.4_ 

Show that if c < —2, then K c does not meet the imaginary axis. 

(Hint: Problem 3.3(a) and Figure 3.7 should help.) 


Problem 3.4 shows that if c < —2, then the set K c is in at least two pieces, and 
so (as when c > |) it is not connected. We pursue the question of the 
connectedness of K c in the next section. 
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4 THE MANDELBROT SET 


SO! 


After working through this section, you should be able to: 

(a) understand the definition of the Mandelbrot set M ; 

(b) use the Fatou-Julia Theorem and its corollaries to determine whether 
certain points lie in M; 

(c) appreciate how a computer may be used to plot M; 

(d) show that certain points c lie in M because the corresponding function P c 
has an attracting cycle ; 

(e) appreciate where certain periodic regions of the Mandelbrot set are located 
by making use of saddle-node bifurcations and period-multiplying 
bifurcations. 


4.1 What is the Mandelbrot set? 

In Section 2 we gave computer-generated plots of a number of keep sets K c for 
quadratic functions of the form P c (z) — z 2 + c. The shapes of these sets are 
remarkably varied, but they can be classified into two distinct types: those 
which are ‘all in one piece’ (for example, K 0 and K- 2 ) and those which are ‘in 
more than one piece’ (for example, K\). The mathematical name for a set 
which is ‘all in one piece’ is connected. 

Earlier in the course we introduced pathwise connectedness , in order to define Unit A3, Section 4 

the key concept of a region. It is clear that the sets Kq and AT_ 2 are both 

pathwise connected, but it is not so evident that the complicated set K -1 is 

pathwise connected, even though it does appear to be in one piece. We find it 

convenient to introduce the following more general notion of connectedness. 


Throughout the rest of this 
unit, this is the meaning of the 
word connected. 


For example, for each real c > |, the set K c is disconnected because, by 
Theorem 3.1 and Problem 2.3(a), it does not meet the real axis, but it does 
have points in both the upper and lower open half-planes. Thus, for c > the 
definition of disconnected is satisfied with A — K c , G\ = {z:\mz > 0}, and 
Gi = {z : Imt < 0}; see Figure 4.1, where c = 5. 

Problem 4.1- 

Use the result of Problem 3.4 to show that K c is disconnected for real c < —2. 



Definitions A set A is disconnected if there are disjoint open sets G i 
and G2 such that 

A l~l Gi ^ 0, A D C?2 ^ 0 and A C Gi U G 2. 

A set A is connected if it is not disconnected. 


It is more difficult to prove that a connected set is connected. This is because 
we need to show that it is not disconnected, that is, the open sets G\, G% in 
the above definition do not exist. However, it is not too difficult to show that 
any pathwise connected set is connected 

(although we omit the details), and so the sets K 0 and K- 2 are connected. At 
present, we set aside the difficulty of proving that a set is connected and 


Figure 4-1 K 1/2 


However a connected set need 
not be pathwise connected. 
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concentrate" on the set of points c such that K c is connected; this set is the 
celebrated Mandelbrot set. 


Definition The Mandelbrot set is the set M of complex numbers c 
such that K c is connected. 


Remark This definition is often phrased in terms of the connectedness of the 
Julia set J c . By using the fact that K c has no holes in it, it can be shown that 
K c is connected if and only if J c is connected. 

For example, since K 0 and AT_ 2 are connected, we have 0 G M and -2 G M. 
On the other hand, by Problem 4.1 and the discussion preceding it, we have 
c ^ M, for all real O \ and c < —2. 

The definition of M is typical of many in mathematics. It defines the object 
that we are interested in precisely, but as it stands it is difficult to work with. 
For most values of c we have very little idea what K c looks like, let alone 
whether or not it is connected! Fortunately, however, there is a numerical 
method of deciding whether K c is connected, based on the following 
fundamental result, the proof of which is outlined in Subsection 4.3. 


Theorem 4.1 Fatou-Julia Theorem 

For any c G C, 

K c is connected 0 G K c . 

This result states that the keep set K c is connected if and only if the point 0 
does not escape to oo under iteration of P c . For example, 0 lies in both Kq and 
2, which are connected, but 0 does not lie in K\, which is disconnected. 

To see the power of Theorem 4.1, try the following problem, which 
characterizes that part of the Mandelbrot set which lies on the real axis. 

Problem 4.2 __ 

(a) Use Theorem 4.1 and Theorem 3.2 to show that if c G [-2, |], then c G M. 

(b) Deduce that M 0 IR = [—2, |]. 


It is natural to ask why the number 0 appears in this special way in 
Theorem 4.1. The reason, as you will see in Subsection 4.3, is that the 
number 0 is the only critical point of each of the functions P c (z) = z 2 + c; 
that is, it is the only point at which each of the derivatives P' c (z) = 2z vanishes. 
By the Local Mapping Theorem, an analytic function fails to be one-one near a 
critical point, and so such points play a significant role in the function’s 
behaviour. 

By using Theorem 2.3, we can turn the condition 0 G K c in Theorem 4.1 into a 
numerical condition which is easier to check. By Theorem 2.3(d), 

0 eK c => P”(0) G K c , for n = 0,1,2,..., 

and so, by Theorem 2.3(a), 

0 G K c => |P"(0)| < r c , for n = 0,1,2,..., (4.1) 

where r c = \ + J\- f |c|. 


Benoit Mandelbrot was born in 
Warsaw in 1924. His initial 
training as a mathematician 
took place at the Ecole 
Polytechnique, Paris. He started 
working for IBM in 1958 and 
became an IBM Fellow in 1974, 
at the Watson Research 
Institute in New York State. 


Pierre Fatou (1878-1929) was a 
French astronomer, although he 
trained as a mathematician. In 
addition to his work on iteration 
theory, he proved an important 
result about the boundary 
behaviour of complex functions. 
Fatou and Julia (see page 27) 
proved (independently) a more 
general version of this result in 
1918-19. After their pioneering 
work in complex iteration, there 
were only a few other 
developments in this field until 
the explosion of interest in the 
1980s, sparked off by the use of 
computers. 


Unit C2, Theorem 3.2 
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Since P c (0) = c, the right-hand side of Implication (4.1) shows that |c| < r c and 
hence, from Figure 2.2, that r c < 2. Thus 

0 e K c =► |P c n (0)| < 2, for n = 1,2,... . (4.2) 

Now, if c satisfies the inequalities on the right-hand side of Implication (4.2), 
then the sequence {P”(0)} is bounded and so 0 e K c , by the definition of K c . 
Hence, by the definition of M and Theorem 4.1, 

c € M -<=> K c is connected 
«=> 0€ K c 

4=» |P"(0)| < 2, for n = 1,2,..., 

and so we obtain the following corollary to Theorem 4.1. 


Corollary 1 The Mandelbrot set M can be specified as follows: 
M = {c : |P”(0)| < 2, for n = 1,2,...} . 


Remark Note that if c $ M, then 0 ^ K c , by Theorem 4.1, and so 
P c n (0) —♦ oo as n —» oo. 

Corollary 1 provides a numerical criterion for determining whether or not a 
point c belongs to M. For any given c, we simply compute the terms of the 
sequence (P”(0)}, which axe 

c,c 2 +c, (c 2 +c ) 2 + c,..., 

and try to decide whether all these terms lie in {z : \z\ < 2 }. 

For example, if c = 0, then {P”(0)} is the constant sequence 

0 , 0 , 0 ,...; 

since all these terms lie in { z : \z\ < 2}, we deduce by Corollary 1 that 0 € M. 
On the other hand, if c = 1, then the terms of {P™(0)} are 

1,2,5,26,...; 

since these terms do not all lie in {z : \z\ < 2}, we deduce by Corollary 1 that 
1 i M. 

Problem 4.3_ 

Use Corollary 1 to determine which of the following points lie in M. 

(a) c = —2 (b) c=l-H (c) c = i (d) c=y/2i 


Corollary 1 makes it possible to use a computer to plot an approximation to 
M. A naive algorithm involves checking the inequality 

|P C "(0)| < 2 (4.3) 

for a large number of points c, and for n = 1,2,..., N, where AT is a suitably 
large positive integer. If Inequality (4.3) is false for some n, then the 
corresponding c lies outside M, but if it is true for n = 1,2,..., N, then c must 
be in M or ‘close to M\ The set M was first plotted in 1979 using an algorithm 
of this kind by Mandelbrot (who had previously been plotting Julia sets). 


For n = 1,2,..., 

Pc n+1 (0) = (P c n (0)) 2 +c. 


A plot of the ‘periodic regions’ 
of M (see Subsection 4.2) was 
made in 1978 by R. Brooks and 
J. P. Matelski. They had 
encountered the iteration of 
quadratic functions while 
studying various groups of 
Mobius transformations. 


37 






A rendering of the set M using this approach is shown in Figure 4.2. As you 
can see, the set appears to be very complicated, consisting of many ‘blobs’ (the 
main one of which is bounded by a cardioid), which Mandelbrot called atoms, 
arranged in a highly organized manner. Some of these atoms are stuck together 
to form a complex molecule, whereas others appear to float free. 



We can use Corollary 1 obtain a number of basic properties of M. 


Corollary 2 The Mandelbrot set M 

(a) is a compact subset of {c: |c| < 2}; 

(b) is symmetric under reflection in the real axis; 

(c) meets the real axis in the interval [-2, |]; 

(d) has no holes in it. 


Proof First note that each term of the sequence {P”(0)} defines a polynomial 
function of c. Indeed 

P c (0) = c, 

P 2 (0) = c 2 +c, 

P 3 (0) = (c 2 + c) 2 + c = c 4 + 2c 3 + c 2 + c, 
and, in general, P"(0) takes the form 

P”(0) = c 2 " 1 + 2 n-2 c 2n 1-1 H-l-c 2 + c, for n > 1. (4.4) 

(a) To prove part (a), we define the sets 

M n = {c: |P"(0)| < 2}, for n = 1,2,..., 

so that Mi = {c: |c| < 2}, M 2 = {c : |c 2 + c| < 2}, and so on. Then, by 
Corollary 1, 


and, in particular, M C M\ = {c:\c\ < 2}. Thus M is bounded. 

Each of the sets M n is closed, because its complement 
C - M n = {c: |P"(0)| > 2} 

is open. Indeed, if |P™(0)| > 2 for some Co, then this inequality must hold 
for all c in some open disc with centre Co, by the continuity of the function 
c i—> |P”(0)|. It follows that M itself must be closed, because if c £ M, 
then c £ M n for some n, and so some open disc with centre c must he 
outside M n and hence outside M. Hence M is closed and bounded; that is, 
M is compact. 


The box represented in 
Figure 4.2 is 

{c : —2 < Rec < 1, 
-1.5 <Imc< 1.5}. 


Note that M is in the c-plane, 
often called the parameter 
plane, whereas each keep set K c 
lies in the 2-plane, often called 
the dynamical plane. 


Statement (d) means that 
C — M is connected. 


This proof may be omitted on a 
first reading. 


This form of P“(0) can be 
justified by Mathematical 
Induction. 


See the proof of Theorem 2.3(b) 
for similar reasoning. 
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(b) Because P"(0) is a polynomial in c with real coefficients, 

Pc"(0)=I>Ro) =S> |i»(0)| = |P c n (0)|, for n = 1,2,... . 

Hence, by Corollary 1, c € M if and only if c G M, and so M is symmetric 
under reflection in the real axis. 

(c) You proved this in Problem 4.2(b). 

(d) The proof that C - M is connected is very similar to the proof of part (f) 
of Theorem 2.3, using {c: |c| > 2} instead of {z : |z| > r c ] and applying the 
Maximum Principle to the analytic function c i—> P"(0) instead of to 

z i —* P?(z)- We omit the details. ■ 

It turns out that the picture of the Mandelbrot set in Figure 4.2 is misleading. 
Some parts of M are so thin that the naive algorithm fails to detect them. This 
became clear when A. Douady and J. H. Hubbard proved the following 
remarkable result in 1982. 


Theorem 4.2 The Mandelbrot set is connected. 


Thus it follows that all parts of M in Figure 4.2 are actually linked together. 
Knowing this fact, we can devise more effective algorithms for plotting better 
renderings of M, such as the one shown (with various enlargements 
corresponding to the small square boxes) in Figure 4.3. 



(a) (b) 


Figure \.3 M — showing its connectedness 

To gain some insight into why the connectedness of M is so remarkable, we 
note that the sets M n = {c: |P"(0)| < 2}, n = 1,2,..., defined in the proof of 
Corollary 2, are in fact nested, that is, 

Mi D M 2 D Mj D • • ■ 

We can represent these nested sets M n readily using a computer. See 
Figure 4.4, where, for the first few values of n, points of M n — M n+ i are plotted 
black if n is odd and white if n is even. Since c € M n — M n+ i if and only if 

|P c fc (0)| <2, for k = 1,2,.. . ,n, but |P c n+1 (0)| > 2, 

these bands measure how long the corresponding sequences {P”(0)} remain in 
{z : \z\ < 2}. It appears that each of the boundaries 

dM n = {c: |P c n (0)| = 2}, n = 1,2,..., 

forms a simple-closed smooth path (that is, it does not break up into several 
pieces). This surprising fact can be shown to be equivalent to the 
connectedness of M. 


It was Douady and Hubbard 
who named the set M after 
Mandelbrot. 


We make no attempt to prove 
Theorem 4.2. 


We discuss the structure of M 
in Subsection 4.2. 



Figure 4.4 


By Equation (4.4), 
aM 1 = {c:|c| = 2}, 

8M 2 = {c:\c 2 + c\ = 2}, 
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4.2 Inside the Mandelbrot set 

Corollary 1 of Theorem 4.1 provides a good means of showing that a given 
point c is not in the set M. This corollary is not so helpful, however, as a 
means of checking that c is in M. Of course, for some values of c, such as c = 0 
or c = -2, it is possible to check directly that |P"(0)| < 2, for n = 1,2,..., but 
this is usually not the case. Instead, the following general result can often be 
used. 


Theorem 4.3 If the function P c has an attracting cycle, then c € M. 


Remark One way to prove Theorem 4.3 is to show that if P c has an We outline another proof of 

attracting cycle, then the sequence {P”(0)} is attracted to this cycle in the Theorem 4.3 in Subsection 4.3. 
way described in Remark 1 following Theorem 2.5. Since {P c n (0)} can be 
attracted in this way to at most one cycle, it follows that P c has at most one 
attracting cycle for each value of c. 

Theorem 4.3 allows us to identify various key parts of M, by determining for 
which values of c the function P c has an attracting p-cycle for various values of 
p. The solution for attracting fixed points and attracting 2-cycles is 
particularly elegant, as we now show. 


Theorem 4.4 

(a) The function P c has an attracting fixed point if and only if c satisfies 

(8|c| 2 — |) 2 + 8Rec < 3. (4.5) 

(b) The function P c has an attracting 2-cycle if and only if c satisfies 

|c+l|<i. (4.6) 


As you will see in the proof of Theorem 4.4(a), the Condition (4.5) is 
equivalent to the statement that c lies inside the cardioid with paxametrization 

7« = - ~e 2it (t e [—7r, 7r]). 

This cardioid is the boundary of the main ‘atom’ of M. Condition (4.6) means 
that c lies inside the open disc {c:|c+l|<|} which lies immediately to the 
left of the cardioid; see Figure 4.5. Thus if a point c lies in one of these two 
sets, then it is possible to verify this using Theorem 4.4, and hence to show 
that c 6 M. For example, the point c = shown in Figure 4.5, seems to 

lie just inside the cardioid. For this value of c, \c\ 2 = | and Rec = — and so 

(8|c| 2 -§) 2 +8Rec=(4-§) 2 —8-f 
= f <3. 

Hence P c has an attracting fixed point, by Theorem 4.4(a), and so c = — | + 
lies in M, by Theorem 4.3. 



An alternative approach is to 
show that |Pc(a)| < 1 for one of 
the fixed points a of P c , but it 
is messier. 


Problem 4.4_ 

Prove that each of the following points lies in M. 
(a) c=-0.9 + 0.If (b) c = 0.2 + 0.5i 
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Proof of Theorem 4.4(a) 

First note that a is a fixed point of P c if and only if 
P c (q) = a 2 + c = q, 
that is, if and only if 
c = a — a 2 . 

Moreover, this fixed point is attracting if and only if 

|P'(a)| = |2 q| < 1. 

Thus P c has an attracting fixed point if and only if c is of the form a — a 2 , 
where |a| < i, that is, if and only if c lies in the image of the open disc 
{z ■ \z\ < 5} under the function z 1—> z - z 2 . To understand the nature of this 
image we use the approach of Unit Dl, Section 4, expressing the function 
z 1 —» z - z 2 as a composition of one-one conformal mappings. In fact, 

and so the function z 1 —► z - z 2 has the effect indicated in Figure 4.6. 



Figure 4.6 


To complete the proof of part (a), we need to show that c lies in the region <S, 
bounded by the cardioid if and only if 

( 8 |c | 2 -§) 2 + 8Rec<3. 

This is a rather fiddly technical exercise, which can be done in various ways. 
One approach is to note that the cardioid itself is the image of the circle 
|z| = | under the function z 1 —> z - z 2 , that is, it is the path in the c-plane 
with parametrization 

7W = \* il ~ (I®") 2 

= \{cost + isint) - |(cos2t + isin2t) (t € [- 7 r, 7 r]). 

In Unit A2, Exercise 2.6, we found that this path has equation 
4|c | 4 — §|c | 2 + 5 Rec = gj, 
which can be rearranged to give 
( 8 |c | 2 — |) 2 + 8 Rec = 3. 

It can then be shown that, for \ < r < |, the cardioid splits the part of the 
circle |c| = r in the open upper half-plane into two arcs (see Figure 4.7). On 
the left arc, which is inside the cardioid, Inequality (4.5) holds, whereas on the 
right arc we have the opposite inequality. For other positive values of r, the 
cardioid does not meet this semi-circle, and the direction of the inequality can 
be found by considering its nature at c = ±r. We omit the details. ■ 


We ask you to prove part (b) in 
Problem 4.5. 


Note that the cusp of the 
cardioid arises because the 
function z 1—* z — z 2 is two-one 
near the point \ (see Unit C2, 
Problem 3.3). In particular, this 
function doubles angles between 
smooth paths emerging from |. 



Figure 4 .7 
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Problem 4.5 _ 

(a) Prove that 

P 2 (z)-z = (P c (z)-z)(z 2 + z + c + 1). 

(b) Deduce from part (a) that if c ^ -f, then P c has the 2-cycle Qi,a 2 where 



with multiplier 

(P c 2 )'(a 1 ) = 4a 1 a 2 . 

What happens if c = -|? 

(c) Deduce from part (b) that P c has an attracting 2-cycle if and only if 
|c+l| < i. 


Theorem 4.4 suggests that each atom in the Mandelbrot set has an associated 
period p such that, for each c in the interior of the atom, the function P c has 
an attracting p-cycle. Numerical experiments appear to confirm this and so we 
introduce the notion of a periodic region. 


Definition A periodic region is a maximal region 72 such that, for 
some positive integer p, 

the function P c has an attracting p-cycle, for all c (E 72. (*) 


For example, the inside of the cardioid and the open disc specified by 
Inequalities (4.5) and (4.6) are both periodic regions. Unfortunately, none of 
the other periodic regions in M seems to have such a straightforward 
characterization, but we can obtain some information about their location by 
using the following result. 


Theorem 4.5 The function P c has a super-attracting p-cycle if and only 
if 

PP( 0) = 0, butP c fe (0)^0, for Jfe = 1,2,... ,p — 1. (4.7) 


Proof By Theorem 2.4, any p-cycle of P c 
a,P c (a),...,P c p - 1 (a), 
has multiplier 

P'(a)P'(P c (a)) • ■ P^Pr 1 ^)) = (2a)(2P c (a)) • ■ • (2Pr 1 (a)), 

since P' c {z) = 2z. Therefore, such a p-cycle is super-attracting if and only if one 
of the points of the p-cycle is 0. But P c has a p-cycle including 0 if and only if 
Condition (4.7) holds, and so the proof is complete. ■ 

For example, if p = 1, then Condition (4.7) becomes 
P c (0) = c = 0, 

as expected, since Po(z) = z 2 has the super-attracting fixed point 0. 

If p = 2, then Condition (4.7) becomes 

P 2 (0) = c 2 + c = 0, but P c (0) = c^0. 

The only solution is c = —1, as expected, since P-i(z) = z 2 — 1 has the 
super-attracting 2-cycle 0,-1. 


This proves Theorem 4.4(b). 


The word ‘maximal’ signifies 
that there is no region satisfying 
(*) which contains 72. but which 
is not equal to 72. 


42 




Problem 4.6_ 

Show that the function P c has a super-attracting 3-cycle for precisely three 
different points c, one of which lies in the interval [—1.8, —1.7] and the other 
two of which form a pair of complex conjugates. 


It can be proved that if P c has an attracting p-cycle, then c lies in a periodic 
region which contains exactly one point Co, say, with a super-attracting p-cycle. 
We call Co the centre of the associated periodic region. Figure 4.8 shows the 
approximate location in the Mandelbrot set of all points Co for which P Co has a 
r super-attracting p-cycle, for p — 1,2, 3,4. 



Figure 4-8 


In fact the boundaries of all the periodic regions are either (roughly) circular in 
shape or are shaped like a cardioid but they are linked together in a very 
complicated manner. We cannot hope to give a full explanation of the way in 
which the periodic regions of M fit together, but we can gain some insight into 
their structure by looking more closely at the two sets given by Theorem 4.4. 
These are plotted in Figure 4.9. 

The key points to consider here are c = which is the cusp of the cardioid, 
and c = — |, where the cardioid and the circle meet. These are the points 
where the cycle structure of P c changes; for example, as c passes through — | 
from the cardioid into the disc, the attracting fixed point of P c becomes 
repelling and the repelling 2-cycle of P c becomes attracting. At such points, we 
say that a bifurcation occurs. In order to characterize such bifurcations, we 
look more closely at the cycle structure of P c for c = | and c = — |. 

For c = j, the function P c = P\/ 4 has just one fixed point, a = |, with 
multiplier 

P» = P{/ 4 (|)= 2 (*) = 1 . 
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For c = — |, the function P c = P-z/a has two fixed points, one of which is at 
a = — |, with multiplier 

^) = ^3/4(-|)=2(-|)=-l. 

Thus, for both these key points c, the function P c has a fixed point a whose 
multiplier P' c (oi) is a root of unity. In order to state a general result, we say 
that the number A is a primitive nth root of unity if A is a root of uni ty 
and if n is the smallest positive integer for which A n = 1. For example, -1 is a 
primitive square root of unity, but it is not a primitive fourth root of unity. 

The following description gives the two types of bifurcation which occur when 
the multiplier of a cycle is a root of unity. 


Theorem 4.6 Suppose that the function P Co ,cq € C, has a p-cycle whose 
multiplier A is a root of unity. 

(a) Saddle-node bifurcation at c 0 If A = 1, then c 0 is the cusp of a 
cardioid-shaped periodic region 72, such that 

P c has an attracting p-cycle, for c € 72. 

(b) Period-multiplying bifurcation at c 0 If A is a primitive nth root 
of unity, for n > 1, then there are two periodic regions 72i and 72 2 
whose boundaries meet at cq such that 


P c has an attracting 


p-cycle, 

np-cycle, 


for c € 72i, 
for c € 72 2 . 


In the next two problems, which are quite challenging, we ask you to 
investigate specific examples of these two types of bifurcations. 


Problem 4.7_ 

(a) Use your solution to Problem 2.2(a) to determine a polynomial Q c (z) such 
that 

P 3 c (z)-z = (P c (z)-z)Q c (z). 

(b) Verify that 

Q_ 7/4 (z) = (z 3 + Iz 2 -fz-|) 2 . 

(c) Deduce that a saddle-node bifurcation occurs at c— — ?, and relate this 
fact to the solution of Problem 4.6 and to Figure 4.8. 


Problem 4.8_ 

Show that if c = C - C 2 , where 2£ is a root of unity (^1), then a 
period-multiplying bifurcation occurs at c. Relate this fact to Figure 4.8, with 
C= -§> ie 2 "/ 3 and \i. 

(Hint: First check that £ is a fixed point of P c .) 


We are now in a better position to describe the structure of the Mandelbrot 
set M. Using the approach of Problem 4.8, we find that each part of the main 
cardioid in Figure 4.8 is decorated by periodic regions. In a similar way all 
these periodic regions are themselves decorated everywhere by further periodic 
regions, and so on. 

All these period-multiplying bifurcations go a long way towards explaining the 
complicated structure of M. In addition, however, we find throughout the 
boundary of M the appearance of small cardioid-shaped periodic regions, 
arising from saddle-node bifurcations, such as the one at c = — | in 
Problem 4.7. All these cardioid-shaped regions are themselves decorated with 
smaller periodic regions, as a result of period-multiplying bifurcations, and so 
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— 1 is the multiplier we found 
above for c = -1 and a = — ±. 


We omit the proof of this 
theorem. 

The name saddle-node 
bifurcation arises from the 
shape of the graph when such a 
bifurcation occurs within a 
family of real functions. 

If n = 2, the name 
‘period-doubling bifurcation’ is 
also used. 


The rest of this subsection is 
intended for reading only. 



they give rise to small copies of the Mandelbrot set within itself (see 
Figure 4.3(c)). These copies appear to be linked together in M (remember that 
the Mandelbrot set is connected) by a complicated network of ‘veins’. The 
simplest such vein lies along the real axis from —2 to A. In fact, one of the 
unsolved problems about the Mandelbrot set M is to decide whether M is not 
only connected, but also pathwise connected. 


The real bifurcation diagram 

Finally, we describe a very effective graphical method of obtaining information 
about the periodic regions of the Mandelbrot set whose centres lie on the real 
axis. The idea is to choose a large number of real values of c between —2 
and |, and plot each of the corresponding sequences {P”(0)} vertically above 
and below a horizontal c-axis (see Figure 4.10(a)). In order to determine any 
attracting p-cycles (to which this sequence will be attracted, by the remark 
following Theorem 4.3) the first 200 or so terms are discarded and the next 600 
or so are plotted. Thus if P c has an attracting p-cycle, then p points should be 
plotted above and below the corresponding point c. 


A pathwise connected set must 
be connected, but a connected 
set need not be pathwise 
connected. However, for open 
sets the two types of 
connectedness are equivalent. 



(a) -2 < c < i (b) -2 < c < -1.54 


Figure 4-10 

As expected, this ‘bifurcation diagram’ reveals the convergence of {P"(0)} to 
an attracting fixed point for -f < c < A, to an attracting 2-cycle for 
— |<c<— |, to an attracting 4-cycle for c just to the left of — |, with further 
period-doubling bifurcations to the left of this. Also visible is an attracting 
3-cycle for c just to the left of — For other values of c, it is less clear what is 
happening, but by scaling the c-axis appropriately (see Figure 4.10(b)) many 
other ‘periodic windows’ are revealed which correspond to attracting p-cycles, 
some of which are labelled in the diagram. Notice that in each such window 
there is a point c whose cycle includes 0; this value of c is the centre of the 
corresponding periodic region in the Mandelbrot set. In fact, pictures of this 
type for the related family of iteration sequences given by 

x n+ i = kx n ( 1 — x n ), with xo = |, where 0 < k < 4, 

were studied in the early 1970s, before the Mandelbrot set itself had been 
plotted. In particular, the order in which the periodic windows appear was 
found and the rate at which periodic-doubling occurs was discovered (by 
M. Feigenbaum) to have a certain universal property. Thus, even the part of 
the Mandelbrot set which lies on the real axis is extremely complicated, and 
there are still unanswered questions about it. For example, it is believed that 
every non-empty open interval of [—2,1/4] meets at least one of the periodic 
windows, but at the time of writing (1993) this has still not been proved. 
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4.3 Outline proofs of Theorems 4.1 and 4.3 

The aim of this subsection is to indicate, without going into all the details, why This subsection may be omitted 
Theorems 4.1 and 4.3 are true. on a first reading. 


Theorem 4.1 Fatou-Julia Theorem 

For any c € C, 

K c is connected «=>■ 0 € K c . 


Theorem 4.3 If the function P c has an attracting cycle, then c € M. 


We shall need the concept of the preimage set P c 1 (E) of a set E under the 
function P c . This is simply the set of points which are mapped to E by P c : 

Pc HE) = {Z : Pc(z) € E}. 

For example, if E = {-1} and c = 0, then PJ' 1 ({— 1 }) = {i, —i}. Note that 
P~ 1 (E) is always symmetric under rotation by n about 0, since P c is an even 
function. 

We say that a compact disc is a compact set whose boundary is a 
simple-closed smooth path (see Figure 4.11). Note that a compact disc need 
not be a disc, but it must be connected. The nature of the preimage set of a 
compact disc E under P c is given by the following result (see Figure 4.12). 


This notation does not imply 
that P c has an inverse function. 



Figure 4-11 


Lemma 4.1 If E is a compact disc and c ^ dE, then P c 1 (E) is 

(a) one compact disc containing 0 , if c 6 int E ; 

(b) two compact discs, neither containing 0, if c 6 ext E. 


The interior and exterior of a 
subset of C are defined in 
Unit A3, Section 5. 



Figure 4-1% 


To indicate why Lemma 4.1 is true we have shown in both parts of Figure 4.12 
a typical ray emerging from the point c and meeting the set E, as well as the 
preimage of this ray which consists of the two rays emerging from 0, combining 
to form a line through 0. 

Now we start the outline proof of Theorem 4.1 by choosing a closed disc 
Eq = {z : \z\ < r}, where 

r > max{r c , |c|} (4.8) r c = I + y/\ + |c| 

and 

Pci 0) $ dE 0 , for n = 1,2,... . (4.9) 

Then we define the sequence of successive preimage sets of Eq under P c : 

Ei = P c 1 (E 0 ), i?2 = P~ 1 (Ei), and so on, Note that the sets Eq,E\,Ei, ... 

are not escape sets here. 
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that is, we put 

E n = {z : P"(z) e E 0 } , for n = 1 , 2 ,... . 

If P c n+ 1 (z) € E 0 , then we must have P"(z) € E 0 also (by Theorem 2 . 2 , because 
r > r c ). Hence 

E n+1 CE n , for n = 0,1,2,..., 
so that the sets E n are nested. Moreover, 

E 0 n Pi n E 2 n ■ • ■ ={z: P c n (z) e E 0 , for n = 0,1,2, .. .} 

= {z : |P c n (z)| < r, for n = 0, 1 ,2,...} 

= K c , (4.10) 

by Theorem 2.2. Thus the shape of K c is determined by the shapes of the 
sets E n . 

Figure 4.13(a) and (b) illustrate the first few of these nested sets E n , 
n = 1,2,..., for the two cases c = —1 and c = 1 . In both cases we have taken 
Eq — {z : \z\ < r}, where r = 1 . 8 , so that Conditions (4.8) and (4.9) are 
satisfied in each case. Points of E n - E n+l are plotted black if n is even and 
white if n is odd. 



(a) K c :c = —1 (b) K c :c= 1 

Figure 4-13 


In Figure 4.13(a), all the sets E n are compact discs. However, in Figure 4.13(b) 
the sets Eo and E\ are compact discs, but E2 consists of two compact discs, P 3 
consists of four compact discs, and so on. We now show that these different 
structures are related to whether or not the point 0 (and hence c) lies in K c . 

First assume that 0 € K c , so that c = P c (0) € K c . Then 


c€intP„, for n = 0 , 1 , 2 ,..., 
so that, by Lemma 4.1(a), 

E n+ 1 = P~ 1 (E n ) is one compact disc, for n = 0,1,2,... . 

Thus E n , n = 0,1,2,..., is a nested sequence of connected compact sets (as in 
Figure 4.13(a)), and it follows from this (though we do not give the details) 
that K c = Eo fl Ei fl • • • is also connected, as required. 

Next assume that 0 £ K c , so that c = P c (0) £ K c . Since c € int E 0 , there is a 
positive integer m such that 

c € int E m - 1 , but c ^ E m . 

Thus, by Lemma 4.1(b), E m is one compact disc, but E m+ 1 is two; see 
Figure 4.14 (which is based on Figure 4.13(b)), in which c = 1 and m = 1 . 


Note that 

c i dE n , for n = 0,1,2,..., 
since P c n (c) = P" +1 (0) £ dE 0 , 
for n = 0,1,2,..., by 
Condition (4.9). 



Figure 4.14 
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Thus, the set K c lies in the union of the two halves of E m+ 1 , and it has points 
in each half because both K c and E m+ i are symmetric under rotation by 7 r 
about 0. Hence K c is disconnected, as required. This completes the outline 
proof of Theorem 4.1. 

Notice from Figure 4.14 that c lies outside both halves of E m+1 and so, by 
Lemma 4.1(b), the set E m+ 2 = P~ 1 (E m+ 1 ) consists of four compact discs (two 
for each half of E m+ i). In general, E m+n consists of 2” compact discs. Now, 
K c is a subset of each of these preimage sets E m+n , and so any connected 
subset K c of K c must lie in exactly one of the 2” compact discs comprising 
E m+n , say E m+n . Thus 

E m +1 2 E m+2 D -DK C . 


It can be proved that the diameters of these compact discs E m+n tend to zero 
as n —► oo, so that K c must be a singleton set, and there will be one such 
singleton set for each of the (infinitely many) possible nested sequences of 
compact discs E m+n . A set obtained by a construction of this type is called a 
Cantor set, and so K c is such a Cantor set for each c not in M. 

It follows that if c ^ M, then K c has no interior points and so the function P c 
cannot have an attracting cycle, in view of Theorem 2.5(a). This shows one 
possible way to prove Theorem 4.3. 


This requires a rather tricky 
application of the Riemann 
Mapping Theorem and 
Schwarz’s Lemma. 

Georg Cantor (1845-1918) 
developed the foundations of 
‘set theory’. For example, he 
showed that the set of rational 
numbers is countable (that is, 
its elements can be arranged to 
form a sequence), but the set of 
real numbers is uncountable. 
Any Cantor set is uncountable. 
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5 BEYOND THE MANDELBROT SET 


After reading through this section, you should be able to: 
(a) appreciate the universal nature of the Mandelbrot set. 


In Subsection 1.4 we discussed briefly the Newton-Raphson function 


N(z) = 


2z 3 + l 


(5.1) 


corresponding to the cubic polynomial function p(z) = z 3 — 1. 


Under iteration of N , all points of C are attracted to one of the zeros of p, or 
else they remain on the common basin boundary (see Figure 1.9), which 
includes the point at oo for the extended function N. 


Following the discovery of the Mandelbrot set, the Newton-Raphson method 
for a general cubic function was investigated by computer, in the early 1980s. 
For most cubic functions, the corresponding Newton-Raphson function behaves 
under iteration in the same way as the function in Equation (5.1), but for some 
cubic functions a difference was found. 


To make this difference precise, we consider the family of cubic functions given 
by 


p c (z) = {z- 1 )(z + \ ~ c)(z + | + c) 

= z 3 -(§ +c 2 )z-!+c 2 , 

where c € C. The corresponding Newton-Raphson function is 

jVc (4 =z -P£(f) = 2z3+ (i- c2 ) 

CK> £(*) 3*»-(!+<?)’ 


(5.2) 


and it is a straightforward matter to check that the critical points of N c (that 
is, the points where AT' vanishes) are the three zeros of p c and the point 0. 
Usually, the critical point 0 is attracted to one of the zeros of p c under iteration 
of N e , although it may also remain on the basin boundary (for example, if 
c = ±(\/3/2)i, then N c ( 0) = oo and N c ( oo) = oo). For some values of c, 
however, the function N c has an attracting p-cycle, where p > 1, to which the 
point 0 is attracted. In Figure 5.1(a), we have plotted in the parameter plane 
those values of c for which the sequence {1V”(0)} does not converge to one of 
the zeros of p c . 




If we zoom in on various parts of this set, then we find small copies of the 
Mandelbrot set! Inspired by examples of this type, Douady and Hubbard 
developed a theory which shows that copies of the Mandelbrot set appear in 
the parameter plane whenever we consider the iteration of suitable families of 
analytic functions. Thus the Mandelbrot set has a universal nature! 


This section is intended for 
reading only. 


For example, if c = ±(\/3/2)i, 
then 

Pc(z) = z 3 - 1. 



(c) 


The boxes represented in 
Figures 5.1(a) and 5.2(a) are 

{c: —2 < Rec,Imc < 2}. 


49 






To emphasize this universal nature, we describe one relation of the Mandelbrot 
set, which is obtained by iterating non-analytic functions of the form 

f c (z) = z 2 +c, 

where c € C. By analogy with Theorem 4.1, Corollary 1, we plot the set 
{c:|/ c n (0)| <2, for n — 0,1,2,...}, 
which is called the tricorn, or Mandelbar set; see Figure 5.2(a). 




This set is symmetric under rotation by 2n/3 about 0 (as well as under 
reflection in the real axis), and it appears to be connected. Closer inspection 
reveals that some parts of the boundary of the tricorn are smooth (for example 
near c = |), whereas other parts are extremely irregular (see Figure 5.2(b)). As 
you might expect, the tricorn contains small copies of itself, but it also contains 
small copies of the Mandelbrot set — a truly universal object! 

In this unit we have only scratched the surface of the subject of complex 
iteration. For example, there is much more to be said about the structure of 
the individual Julia sets J c , and there are many families of entire functions 
(such as z i—> e cz , where c is a complex parameter), whose behaviour under 
iteration leads to completely new phenomena. Nevertheless, we hope that you 
have gained some insight into this remarkable subject, and that you appreciate 
the irony in the following quotation from A. Douady. 

I must say that, in 1980, whenever I told my friends that I was just 
starting with J. H. Hubbard a study of polynomials of degree 2 in one 
complex variable (and more specifically those of the form z i—> z 2 + c), 
they would all stare at me and ask: Do you expect to find anything new? 



(c) 


From H. -O. Pietgen and P. H. 
Richter, The Beauty of Fractals 
(Springer-Verlag, 1986). 
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EXERCISES 


Section 1 

Exercise 1.1 Plot the terms zo, Zi, 2 2 > 23 for each of the following iteration 
sequences, and write down the corresponding function /. 

(a) z n+ 1 = iz n , z 0 = 2i 

(b) z n+ 1 = 2z n (l - z n ), z 0 = \ 

( C ) ^n+l = ^-j-, 20 = |( —1 + l) 

Exercise 1.2 Find all the fixed points of each of the following functions and 
classify them as (super-)attracting, repelling or indifferent. 

(a) f{z) = z-z 2 (b) f(z) = 2z(l-z) 

(c) /(*) = **-! (d) /(*)=*/(*+ 1) 

Exercise 1.3 Prove that any iteration sequence of the form 
z n+1 = 2z n (l - z n ), n = 0,1,2,..., 
is conjugate, via the conjugating function h(z) = 1 — 2z, to one of the form 
w n+ i=wl n = 0,1,2,.... 

Hence obtain a formula for z n in terms of zo- 

Exercise 1.4 Let f{z) = (az + b)/(cz -I- d), where ad-bc^Q and c ^ 0, let / 
be the corresponding extended Mobius transformation, and let z n = f n {z 0 ), 
where zo € C. 

(a) Prove that if (a — d) 2 + 46c = 0, then / has a unique fixed point a, which 
lies in C. Deduce that the sequence {z n } is conjugate, via the conjugating 
function h(z) = l/(z — a), to an iteration sequence of the form 

w n+1 = w n + 2c/{a + d), n = 0,1,2,.... 

Hence prove that z n —► a as n —► 00 , for all zq € C. 

(b) Prove that if (a — d) 2 + 46c / 0, then / has two fixed points a and /?, both 
lying in €. Deduce that {z n } is conjugate, via the conjugating function 
h(z) = (z — a)/(z — /?), to an iteration sequence of the form 

tu„+i = f'(a)w n , n = 0,1,2,.... 

Hence prove that if |/'(a:)| < 1, then z n —> a as n —> 00 , for all 
z 0 € C — {/?}. Also describe the behaviour of {z n } if 
(i) \f'(a)\ > 1, (ii) |/'(a)| = 1. 

{Hint: You will need to use the fact that a and j3 are fixed points of / to 
simplify the algebra.) 


This exercise, which is fairly 
challenging, contains a complete 
analysis of the behaviour of 
iteration sequences defined by 
Mobius transformations. 
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Section 2 


Exercise 2.1 Use Theorem 2.1 to show that the iteration sequence 
z„+i = 3z„(l - z n ), 71 = 0,1,2,..., 
where zq = |, is conjugate to an iteration sequence of the form 
w n +i — w 2 4- of, n = 0,1,2,... . 


Exercise 2.2 

(a) Prove that if |c| < |, then 



\Pc(*)\ < | + \ 


(b) Deduce from part (a) that if |c| < |, then 
{z:\z\<l + y/\-\c\)cK c . 


(c) Combine the result of part (b) with Theorem 2.3(a) to show that, when c 
is close to 0, the set K c is approximately equal to the closed unit disc. 


Exercise 2.3 For each of the following functions / and points a, show that a 
is a periodic point of / and decide whether it is (super-)attracting, repelling or 
indifferent. 


(a) f(z) = -z, a = i (b) f(z) = z* - 2, <* = ±(-l + v'S) 

(c) f(z) = z 3 + i, a = 0 (d) f(z) = z 3 , a = e ni/13 


Exercise 2.4 This exercise relates to Properties (2.12) and (2.13) for the 
function Po(z) = z 2 , whose Julia set is the unit circle. 

(a) Determine the repelling periodic points of P 0 and prove that each arc of 
the unit circle contains such a repelling periodic point. 

(b) Determine the backward iterates under Pq of the point 1 and prove that 
each arc of the unit circle contains such a backward iterate. 


Section 3 

Exercise 3.1 Let P i/4 (x) = x 2 + ±. 

Plot y = Pi/4,(x) and y = x on the same diagram and use graphical iteration to 
plot the iteration sequence 

x„+i = Pi /4 (x„), n = 0,1,2,..., 
with xo = 0. Describe the behaviour of the sequence {x n }. 

Exercise 3.2 

(a) Use graphical iteration to show that any iteration sequence of the form 

Xn+1 = r~+I’ ” = 0,1,2,..., 

with x 0 € IR - A, where A = {—1, —5, —5, ■ • •}, converges to the point 0. 
(Note that if x = -1/n, then x/(x + 1) = -1 /(n - 1).) 

(b) Relate your answer to part (a) to Exercise 1.4(a). 


In order to attempt this exercise 
you will need to have read 
Subsection 2.4. 
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Section 4 


Exercise 4.1 Use Corollary 1 to Theorem 4.1 to determine which of the 
following points c lie in M. 

(a) c=—l + 2i (b) c=—1 (c) c = —1 + i 

Exercise 4.2 Show that if |c| = 2 but c ^ —2, then 
|c 2 + c| > 2, 

and deduce that c ^ M. 

(Hint: Use the factorization c 2 + c = c(c + 1).) 

Exercise 4.3 Let M n = {c: |P c n (0)| < 2}, for n = 1,2,... . 

(a) Prove that M n C Mi, for n = 2,3,... . 

(Hint: If |c| > 2, then |c| > r c , by Figure 2.2.) 

(b) Prove that M n+1 C M n , for n = 1,2,... . 

(Hint: If |c| < 2, then 2 > r c , by Figure 2.2.) 

Exercise 4.4 Prove that each of the following points lies in M. 

(a) c= -1.1 -0.lt (b) c = 0.6i 

Exercise 4.5 

(a) Prove that if TZ is a periodic region of M, then &R C M. 

(b) Deduce from part (a) that each of the following points lies in M. 

(i) c= \ + (ii) c= -1 + \i 

(Hint: Locate these points in Figure 4.8) 

Exercise 4.6 Explain why there can be at most six values of c (see Figure 4.8) 
for which P c has a super-attracting 4-cycle. 

Exercise 4.7 Show that a period-doubling bifurcation occurs at c— — |, and 
relate this fact to Figure 4.8. 




SOLUTIONS TO THE PROBLEMS 


Section 1 

1.1 (a) i, -1, 1, 1. 



Here /(z) = z 2 . 

(b) 0, 1, |, l 


1.3 (a) If /(z)-=-z- +-b, then 
f'(z) = z + b, 

f\z) = (z + b) + b = z + 2b, 
f(z) = (z + 2b) + b = z + 3b , 
and, in general, 

/ n (z) = z + nb, for n = 1,2, 
(b) If /(z) = z 3 , then 

f'(*) = * 3 > 

/ 2 (z) = (z 3 ) 3 = Z 9 , 
/3(Z)=(Z 9 ) 3 =Z 22 , 
and, in general, 

/ n (z) = z 3 ", for n = 1,2,... 



Here /(z) = \z + 1. 
(c) 0, -1, 0, -1. 


x 



Here /(z) = z 2 — 1. 
(d) 0, i, -1 + t, -i. 


1.4 (a) If zo = 1, then we have 
z n = 1, for n = 1,2,..., 
so that {z n } is constant, and converges to 1. 

(b) If zo = — i, then the terms of the sequence are 

1,1,1,..., 

so that {z„} is eventually constant, and converges to 1. 

(c) If zo = e 2,rl/3 , then the terms of the sequence are 
e 2W3 e 4*«/3 e 8*i/3 _ g 2iri/3 

so that the sequence cycles endlessly between these two 
values. 

(d) Since z„ = z\ and |l/zo| < 1, we deduce that 



is a null sequence. Hence z n —> oo as n —» oo, by the 
Reciprocal Rule (Unit A3 , Theorem 1.5). 



Here /(z) = z 2 + i. 


1.2 If /(z) = \z + 1, then 
f 2 (z) = /(/(*)) 

-/(** + !) 

= Kf*+ 1 ) + 1 


/ 3 (z) = f 2 (f(z)) 


1.5 (a) The fixed point equation is /(z) = |z + 1 = z, 
and 

Iz+l = z <!=> iz = l, 
so the only fixed point is 2. 

(b) The fixed point equation is /(z) = z 2 — 2 = z, and 

z 2 - 2 = z z 2 — z — 2 = 0, 

so the only fixed points are 2 and —1. 

(c) The fixed point equation is /(z) = z 3 = z, and 

* 3 = * <=> z(z 2 -l)=0, 

so the only fixed points are 0, 1 and —1. 

1.6 (a) Since/'(z) = 2z, we have 

|/'(0)| = 0 < 1, so 0 is an attracting (in fact, a 
super-attracting) fixed point of /; 

|/'(1)| = 2 > 1, so 1 is a repelling fixed point of f. 

(b) Since /'(z) = |, we have 

|/'(2)| = \ < 1, so 2 is an attracting fixed point of /. 

(c) Since /'(z) = 2z, we have 

|/'(2)| = 4 > 1, so 2 is a repelling fixed point of /. 
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1.7 (a) From Example 1.2(a), we know that, for 
n = 1 , 2 ,..., 

/”(*) = ^ z 

-*0asn^oo, for all z € C. 

Hence the basin of attraction of 0 under / is 
i z : / n ( 2 ) —► 0 as n —> oo} = C. 

(b) From Problem 1.3(b), we know that, for n = 1,2,..., 

rw = * 3 ". 

f n (z o) -* 0 as n -too, for |z 0 | < 1, 
but 

f n (z o) -** 0 as n —► oo, for |zo| > 1. 

Thus the basin of attraction of 0 under / is 
{ z : f n ( z ) —+ 0 as n —* oo} = {z : |z| < 1}. 


1.8 First we note that the function h(z) = —z + | is 
one-one on C. If w n = h(z n ) = —z„ + f, then 
z„ = -w n + \, so 

Zn+l = Zn-zl , n = 0,1,2,. . . , 
becomes 

—ut n+ i + | = (—w„ + 5) — (—tt; n + |) 

for n = 0,1,2,..., 


that is, 

w n+ i=tt£ + i, for n = 0,1,2,..., 
which demonstrates that these two sequences are 
conjugate. 

If zo = then wo = 0. 


(i) If |a| < 1, then {a n } is a null sequence, so in this case 

b -* 

z„ -> 1 _ as n -t 00. 

(Note that 6/(1 - a) is the only fixed point of the 
function f(z) = az + 6.) 

(ii) If |a| = 1, a / 1, then {a n } is divergent, by Unit AS, 
Theorem 1.7(b). It follows that {z„} is divergent in this 
case (unless zo = 6/(1 — a), in which case {a„} is 
constant). 

(iii) If |a| > 1, then {a n } tends to infinity, by Unit AS, 
Theorem 1.7(a). It follows that {z n } tends to infinity in 
this case (unless zo = 6/(1 - a)). 


1.10 With N(z) = ^ -- and h(z) = - —we have 

2z + a z — (3 

_ z 2 - 6 - a(2z + a) 

- z 2 - 6- 0(2z + a) 

_ z 2 - 2az - (qq + 6) 
z 2 - 2 0z - (a/3 + 6) 

= y~i7 z +} (bythehint) 

= (z-af 
(z — P) 2 
= (h(z))\ 

as required. 

1.11 If p(z) — (z — a) 2 , then 

and so, by Problem 1.9(b)(i) (with a = 5,6 = |a), all 
points of C are attracted to a under iteration of N. 


1.9 (a) First we note that the function 
h(z) = z + b/(a-l), 
where a / 1, is one-one on C. 

If w n = h(z n ) = Zn + 6/(a - 1), then 
Zn = Wn- b/(a - 1), SO 

z n+1 = az n + 6, 11 = 0,1,2,..., 

becomes 


wn+i - b/(a - 1) = a(w n - b/(a - 1)) + 6, 

for n = 0,1,2,; 

that is, 

w„+i = aw n , for n = 0,1,2,..., (1) 

since -b/(a - 1) = -ab/(a - 1) + 6. 

(b) The iteration sequence (1) has general term 
w n = a n w o (see Example 1.2(a)). Hence 


giving 


that is, 



for n = 0,1,2,..., 


for n = 0,1,2,...; 

for n = 0,1,2,... . 


Section 2 

2.1 (a) By Theorem 2.1, with a = —4, 6 = 4, c = 0, the 
sequence z n +i = —4z 2 + 4z n is conjugate to the sequence 

Wn+l =wl+d, 

where d = ac + 56 — |6 2 = 2 — 4 = —2, using 
h(z) = —4z + 2. In this case, wo = = 0. 

(b) By Theorem 2.1, with a = —2, 6 = 0, c = 1, the 
sequence z n+ i = 1 — 2z 2 is conjugate to the sequence 

where d — ac+ \b- \b 2 = -2, using h(z) = -2z. In this 
case, wo = h(0) = 0. 

2.2 (a) P 2 (z) = P c (P c (z)) 

= {z 2 + c) 2 +c 
= z 4 + 2cz 2 + c 2 + c; 

P c 3 (z) = (z 4 + 2cz 2 + c 2 + c) 2 + c 

= z 8 + 4cz 6 + (6c 2 + 2c) z 4 + (4c 3 + 4c 2 ) z 
+ c 4 + 2c 3 + <? H- c. 
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(b) We have 

P? + \z) = -Pc(Pc n (*)) 

= (P c n (z )) 2 + c. ( 1 ) 

Since P c (z) = z 2 + c, it follows from Equation (1) that 
the degree of P " +1 is twice the degree of P", and hence 
the degree of P" is 2 n . The evenness of P" follows from 
the fact that P c is even, because 

P c (-*) = (-z) 2 -t-c = * 2 +c=P c (z). 


Now the function J is one-one on (u>: |ui| > 1 } and hence 
w n+ i=wl, for n = 0 , 1 , 2 ,... . 

Since |tu 0 | > 1, we have tu n -* oo as n -> oo, and hence 
z n = »n + l/w„ —» oo as n —* oo (by the Reciprocal 
Rule). 

(c) By part (a), we find that no point of L belongs to 
E- 2 , whereas by part (b) all points of C - L belong to 
P- 2 - Hence P_ 2 = C — L, and P _2 = C — P_ 2 = L. 


2.3 (a) The fixed point equation is P c (z) = z 2 + c = z, 
and 

Thus the fixed points are 



Let these fixed points be a and 0, so that a + 0 = 1, and 
Pc(a) = 2a and P' c {0) = 20. 

Now a + 0 = 1 and so \{P' c {a) + P'c{0)) = 1- Unless 
c = j (so that Pc (a) = P' c {0) = 1), this implies that at 
least one of P' c {a), P' c {0 ) lies outside the unit circle, and 
so the corresponding fixed point is repelling. 

(In fact, if a = \ + y/\ - c, then P' c {a) lies outside the 
unit circle.) 

(b) If c = j, then P c = P 1/4 has just one fixed point, 1, 
which is indifferent since K/S) = 1 - 

2.4 (a) ro — | ^^ + 0 = 5 + 5 = 1 ; 

r * = 5 + \/s +1*1 = 5 + \[\ = K 1 + v/ 5); 

r -2 = 5 + \Jk + \~ 2 \ = \ + \/f = 2 ‘ 

(b) The sequence {Pq 1 (1)} has terms 1,1, 1 ,..., and so 
P 0 n (l) -** 00 as n —> 00. This shows that ro = 1 is the 
smallest value of r 0 for which Theorem 2.2 holds with 
c = 0 . 

The sequence {P_ 2 ( 2 )} has terms 2,2,2,..., and so 
P -2 (2) -« 00 as n -too. This shows that r _2 = 2 is the 
smallest value of r _2 for which Theorem 2.2 holds with 
c = - 2 . 


2.6 To prove that E c is completely invariant under P c , 
we note that 

2 € P c P c n (2) —> 00 as n —► 00 

^ P c " +1 (z) - 00 as n -* 00 

Pc n (P c (2)) - 00 as n - 00 
<=► ^c(2)€Pc, 

as required. 


2.7 (a) Since 

P i (-t) = (-i) 2 +t = -l+i 
and 

P < (-l + i) = (-l + i ) 2 + i = -i, 

P?(-i) = -i and P 2 (-l + i) = -i + j ) 
and so — i and —1 + i form a 2-cycle of Pi. Hence both 
these points lie in Ki as do i and 1 — i, by 
Theorem 2.3(e). Another point lying in Ki is 0, because 
Pi(0) = i and i e Ki. 




Clearly none of these points axe fixed points of Ki. 
(b) Since Pq(z) = z 8 , we have to solve the equation 


2.5 (a) If zo 6 L, then zo is real with |zo| < 2, and so 
z\ = Zq 2 

is real with —2 < z\ < 2, because 0 < Zq < 4; hence 
zx € L. On repeating this process, we deduce that 
z n S L, for n = 1 , 2 ,... . 


If zo € C — L, then z\ € C — L. For if zi € L, then 
zo = 21 + 2 € [0,4] => 2 0 €L. 

On repeating this process, we deduce that 
2 „eC —L, forn = 1 , 2 ,... . 

(b) Let w n = J _ 1 ( 2 „), for n = 0,1,2,..., so that 
|w„| > 1 , for n = 0 , 1 , 2 ,... . 

Then z„ = J(w n ), and so, for n = 0 , 1 , 2 ,..., the 
equation z n + 1 = z 2 — 2 gives 
J(Wn + l) = (J(Wn)) 2 - 2 


= K+:f -2 


^ z(z 7 - 1 )= 0 . 

The solutions are 0 and e 2lrk 'k = 0,1,..., 6 (Unit Al, 
Theorem 3.1). 

Of these, the points 0 and 1 (k = 0) are fixed points of 
Po, whereas 

Po(e 2 ^ 7 )=e 4 ^ 7 , 

Po(e-' 7 )=e-' 7 , 

Po(e 8 * i/T ) = e 16 " </7 = e 2 ' i/T , 
and 


Po(e 6 * i/7 ) = e 12lr</7 
P 0 (e 12lr< ^ 7 ) = e 24lr< ^ 7 = e 10lr *^ 7 
Po(ei 0lr */ 7 ) — e 20lri P — e 6lT 'P 
Hence 

g 2*</7 g 4 wi/7 g8TTi/7 an j g 67ri/7 g 12«/7 gl0*i/7 




are both 3-cycles of Po. 
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Using this and the results of Example 2.1, we obtain the 
following diagram. 


Section 3 

3-1 (a) 



(c) Since 

P_5/ 4 (^(-l + V2)) = (i(-l + a/ 2)) 2 - | 

= \(3 — 2a/2) - | 

= i(-l-V2) 

^i(-l + V2), 

Pl s , 4 (i(-l + v^)) =P_ s/4 (i(—1 -a/2)) 

= (i(-l - a/ 2)) 2 - | 

= £(3 + 2a/2) - | 

= i(-l + V2), 

the point |(—1 + a/ 2) is a periodic point, with period 2, 
of P_5/ 4 . (It forms a 2-cycle with j(—1 — \/2).) 

2.8 (a) The numbers —i, — 1 + i form a 2-cycle of Pj, 
with multiplier 

(pf)'(-t) = p'H)P'(-i + i) 

= (—2i)(—2 + 2i) 

= 4 + 4 i, 
by Theorem 2.4. Hence 
|(P?)'(-i)| =4V2>1, 
and so this 2-cycle is repelling. 

(b) The numbers e 2lr,/ ' 7 , e 4l " /7 , e 8lr,/7 form a 3-cycle of 
Po, with multiplier 

(P 0 3 )'(e 2 " /7 ) = P^e 2 ^ 7 ) P^e 4 ^ 7 ) P^(e 8 ' t/7 ) 

= (2e 2 ^ 7 )(2e 4 "/ 7 )(2e 8 ^ 7 ) 

= 8 , 

by Theorem 2.4. Hence 
| |> 1 , 
and so this 3-cycle is repelling. 

(c) The numbers |(-1 + a/ 2), |(-1 - \/2) form a 
2-cycle of P_ 5 / 4 , with multiplier 

(P! s/4 )'(i(-l + V2)) = Pi s/4 (i(-l + V2)) 

= (-1 + a/ 2)(-1 - a/ 2) 

by Theorem 2.4. Hence 

|(p! 5/4 ) , (I(-1 + V2))|^1, 

and so this 2-cycle is indifferent. 



y = —2x + 1 


(b) If xo = 0, then {a; n } tends to infinity and if xq = 1, 
then the sequence is constant. In Problem 1.9(b)(iii) we 
found that *„ +1 = az n + 6, n = 0, 1 , 2 ,..., tends to 
infinity if |a| > 1 unless z 0 = 6/(1 - a). Here we have 
a = —2, 6=1 and 6/(1 — a) = |, and so our answers 
agree with this result. 

3.2 (a) With xo = 0, graphical iteration gives the 
following diagram. 



(b) Since 

x 2 + 1 > x <=)> x 2 - x + 1 > 0 

<=> (*-§)’ +f >0, 

we deduce that 

x 2 + 1 > x, for all x e R. 

Hence the graph y = x 2 + 1 lies strictly above y = x, that 
is, the function f(x) = x 2 + 1 has no real fixed points. 
Also x n +i — x„ + 1 > x n and so the sequence {a:,,} is 
increasing and must tend to infinity (because there are no 
fixed points to prevent this). 

(c) Since x„ = P"(a:o), for n = 1,2,..., we deduce that 
Pi(x o) —* oo as n —► oo, for all x 0 € R, 

and so no point of IR belongs to K\. 
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3.3 (a) If c is real and y is real, then 
p c(iy) = ( iy) 2 + c= —y 2 + c, 

which is real. 

(b) Theorem 3.1 implies that if c > then 
P"(x) -too as n -too, for all x € IR, 
and hence, by part (a), that 

Pc(Pc{iy)) -* oo as n -t oo, for all y € IR. 

Thus 

Pc +1 (iy) —♦ oo as n -* oo, for all y G IR, 
and so no point of the imaginary axis belongs to K c . 

3.4 If c < -2 and y E IR, then 
Pc{iy) = -y 2 + c 



and so P c (iy) lies outside the interval I c (see Figure 3.7). 
Hence 

P?(Pc(iy))-> oo as noo, for all y € IR, 
by (3.3), and so 

P? +1 (iy) —♦ oo as n -* oo, for sill y € IR. 

Thus no point of the imaginary axis belongs to K c . 

Section 4 

4.1 In Problem 3.4 we saw that, if c < —2, then K c does 
not meet the imaginary axis. Since K c has points in both 
Gi = {z : Re 2 > 0} and G2 = {2 : Rez < 0} (for example, 
the fixed points | ± y/\ — c), we deduce that K c is 
disconnected for c < — 2. 

4.2 (a) By Theorem 3.2, K c contains 0 if c € [—2, 
and so, by Theorem 4.1, K c is connected and hence 
c€ M. 

(b) Since K c is disconnected for c > | and c < —2, by 
Problem 4.1 and the discussion preceding it, we have 
c ^ M for c > j and c < —2. Hence, by part (a), 

M nlR= [-2,i]. 

4.3 (a) If c = —2, then the terms of {P"(0)} are 
-2,2,2,... . 

Since all these terms lie in {z : \z\ < 2}, we deduce by 
Corollary 1 that —2 € M. 

(b) If c = 1 + i, then the terms of {P"(0)} are 
1+t, (1 + i) 2 + 1 + t = 1 + 3i, .... 

Since |1 + 3i| = %/l0 > 2, we deduce by Corollary 1 that 

1 + i i M. 

(c) If c = i, then the terms of (P c n (0)} are 
i,-l+ *,-*,-1 + t,... . 

Since all these terms lie in {2 : |^| < 2}, we deduce by 
Corollary 1 that t € M. 

(d) If c = \/2i, then the terms of (P"(0)} are 

V2i,-2 + V2i,... . 

Since |—2 + y/2i\ = y/6 > 2, we deduce by Corollary 1 
that V2i £ M. 


4.4 (a) The point c = —0.9 + O.lt appears to lie in the 
disc |z + 1| < so we use Theorem 4.4(b). Since 

- |c+1| =4_Q.1J-0+tl= 0.1414... < i, 

Pc has an attracting 2-cycle, by Theorem 4.4(b). Thus c 
lies in M, by Theorem 4.3. 

(b) The point c = 0.2 + 0.5i appears to lie inside the 
cardioid, so we use Theorem 4.4(a). Since |c| 2 = 0.29 and 
Re c = 0.2, we have 

(8|c| 2 - §) 2 + 8Rec = (2.32 - 1.5) 2 + 1.6 
= 2.2724 < 3, 

and so P c has an attracting fixed point, by 
Theorem 4.4(a). Thus c lies in M, by Theorem 4.3. 

4.5 (a) Since P 2 (z) = (z 2 + c) 2 + c, we have 
P, 2 (z) - z = z 4 + 2cz 2 - z + c 2 + c. 

(P c (z)-z)(z 2 + z + c+l) =(z 2 -z + c)(z 2 + z + c + 
= z 4 + 2cz 2 - z + c 2 + c, 

as required. 

(b) The 2-cycles of P c are the solutions of P 2 (z) — z = 0 
which are not solutions of P c (z) — z = 0. Hence, by 
part (a), they are the solutions of 
z 2 + z + c + 1 = 0; 
which gives the 2-cycle c*i, «2 where 



Note that c = — | must be excluded, because in that case 
oti = «2 = —5, which is a fixed point of P-3/4. By 
Theorem 2.4, the multiplier of this 2-cycle is 

(P c 2 )'(a 1 ) = P'(a 1 )P'(a 2 ) 

= (2ai)(2a 2 ) 

= 4aiQ2. 

(c) Since 0:102 = c + 1, we deduce that the above 
2-cycle is attracting if and only if 

|(f’c 2 )'(<*i)| = |4aia 2 | = 4|c + l| < 1; 
that is, if and only if |c+ 1| < j, as required. 

4.6 The function P c has a super-attracting 3-cycle if and 
only if 

P c 3 (0) = (c 2 + c) 2 + c = 0, 

but P c (0) = c 0, P 2 (0) = c 2 + c ^ 0. Thus we seek the 
solutions of 

c 4 + 2c 3 + c 2 + c = 0 

which are not c = 0 nor c = — 1. Since c = 0 is a solution 
and c = — 1 is not, we need to solve 
c 3 + 2c 2 + c + 1 = 0. 

Using calculus, we see that the real function 

/(c) = c 3 + 2c 2 + c + 1 has the following graph, so that 

there is only one real zero. 





Since /(—1.8) = -0.152 and /(—1.7) = 0.167, it follows 
that this value of c lies in [—1.8, —1.7]. 

The remaining two solutions form a pair of complex 
conjugates. (In fact the three solutions, correct to three 
decimal places, are —1.755 and —0.123 ± 0.745i. (See 
Figure 2.18(b)!) These solutions may be obtained by 
using Formula (0.6) in the Introduction to Unit All) 


4.7 (a) Since any solution of P c (z ) - z = 0 is also a 
solution of P 3 (z) — z = 0, we expect P c (z) — z to be a 
factor of P 3 (z) — z = 0. By inspection, we find that 
P?(z) - z = (P c (z) - z)Q c (z) 

= (z 2 -z + c)Q c (z), 

Qc(z) = Z 6 + z s + (3c + l)z 4 + (2c + l)z 3 

+ (3c 2 + 3c + l) z 2 + (c 2 + 2c + l) z 
+ c 3 + 2c 2 + c + 1. 

(Begin by finding the coefficient of z 6 and the constant 
term, and then work inwards from both ends, obtaining 
the coefficient of z 3 twice, as a check.) 

(b) On substituting c = we obtain, after some 
arithmetic, 

Q—7/4(2) = Z 6 + Z 5 — — f Z 3 + ||z 2 + JgZ + jrj 



by inspection again. 

(c) The graph y = x 3 + |x 2 - |r — | is as follows. 



z 3 + ±z 2 — §z — i = (z — a)(z - 0 )(z - 7), (1) 

where a, /?, 7 are the real zeros in the figure. Since these 
are the only solutions of P 3 7/4 (z) — z = 0, which are not 
fixed points of P- 7 / 4 , they must form a 3-cycle of P- 7 / 4 ■ 
By Theorem 2.4, the multiplier of this 3-cycle is 
(P 3 )'(a) = P'(a)P'(/?)P'(7) 

= (2a)(2/3)(2 7 ) 

= 8a/37 
= 1, 

since a/?7 = |, by Equation (1). 

Hence, by Theorem 4.6(a), a saddle-node bifurcation 
occurs at c = —7/4, so we expect to see a small 
cardioid-shaped periodic region with cusp at this point. 
This is included in Figure 4.8, and its centre is at the 
point c = —1.755, found in Problem 4.6. 

4.8 Since 

P c (0 = c 2 + C = < (c = C-< 2 ), 
we find that £ is a fixed point of P c . Also, the multiplier is 

Pc(0 = 2C, 

which is a root of unity (# 1). Hence, by Theorem 4.6(b), 
a period-multiplying bifurcation occurs at c. (In fact, 
since the cardioid is the image of the circle |z| = 1 under 
the function /(z) = z — z 2 , the point c lies on the 
cardioid.) 

If C = — 5, then 2Q = —1, which is a primitive square root 
of unity and so a period-doubling bifurcation occurs at 

and this is visible in Figure 4.8. 

If C = f e 2 ™/ 3 , then 2£ = e 2?r!/3 , which is a primitive cube 
root of unity and so a ‘period-trebling’ bifurcation occurs 
at 

c = je 2 **/ 3 — (je 2 ”^ 3 ) 2 = — | + |V3», 
and this is visible in Figure 4.8. 

If C = 5*1 then 2( = i, which is a primitive fourth root of 
unity and so a ‘period-quadrupling’ bifurcation occurs at 

and this is visible in Figure 4.8. 
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SOLUTIONS TO THE EXERCISES 


Section 1 

1.1 (a) z 0 = 2i, zi = -2, z 2 = -2i, z 3 = 2. 



Here f(z) = iz. 

(b) 20 = j,zi = §,2 2 = 55,24 = |f|. 


/ // 


Here f(z) = 2z(l - z ). 

(c) z 0 = i(-l + t),2i = i,z 2 = j(l + i),2 3 = 5(2 + i). 



Here /(z) = ITT 

1.2 (a) We have to solve f(z) = z: 
z-z 2 =z <=> z 2 = 0 
<!=► 2 = 0. 

Since f'(z) = 1-22, we have 
/'( 0 ) = 1 , 

so the fixed point at 0 is indifferent. 

(b) We have to solve f(z) = z: 

2 z(l - 2) = 2 <=*> 2(1 - 2 2) = 0 , 

so the fixed points of / are at 0 and 1. Since 
f'(z) = 2 — 42, we have 

/'(0) = 2 and f' (5) = 0. 

Thus the fixed point at 0 is repelling, whereas that at § 
is super-attracting. 

(c) We have to solve f(z) = 2: 

so the fixed points of / are at 1(1 ± \/3). Since 
f\z) = 22, we have 

|/'(i(l +V3))| =l + >/3 = 2.732... 
and 

|/'(i(l-V3))| = 0.732... . 

Thus the fixed point of / at 5(1 + \/3) is repelling, 
whereas that at j(l - \/3) is attracting. 


(d) We have to solve f(z) = 2: 



so the only fixed point of / is at 0. Since 
/'(*) = l/(* + l) 2 . we have 

/'(0) = 1. 

Thus the fixed point of / at 0 is indifferent. 

1.3 Putting w n = h(z n ) = 1 - 2z n for n = 0,1,2,..., we 

2„ = 1(1 — w n ), for n = 0,1,2,... . 

Hence the given sequence is conjugate to the sequence 
{iu„}, where 

1(1 - Wn+1 ) = 2(1(1 - u, n )) (1 - 1(1 - Wn )) 

= i(l-«n)(l + tn„) 

= ^(1-^ 2 ), 

so that w n +i = w 2 , for n = 0,1,2,..., as required. 

Since w n = wf, for n = 0,1,2,..., by Example 1.2(b), 
we deduce that 

= !(i-«») 

= 5 ( 1 - w 2 ") 

= 1^1 -(1 - 22 0 ) 2 "V for n = 0,1,2,... . 


1.4 To find the fixed points of / we solve f(z) = 
so the fixed points of / are at 


i(a-d± > /(a-d) a +46c). 


(a) If (a — d) 2 + 4 be = 0, then / hcis just one fixed point 
at a = (a — d)/(2c). Since 00 is not a fixed point of / 
(because c / 0), the only fixed point of / is at a. Putting 
Wn = h(z n ) = l/(zn - a), we have z n = a + l/w n , for 
n = 0 , 1 , 2 ,... . 

Hence the sequence {2 n } is conjugate to the sequence 
{io n }, where 

_ a(a + l/w n ) + b 
~ c(a + l/w n ) + d 
_ aa + b +a/wn 
ca + d + c/ w n 
_ (aa -(- b)w n 


f 1/Wn+l = - 


Thus 


(ca + d)w n + c 


1 _ (aa + b)w n -I- a — a(ca + d)w n 

Wn+l (ca + d)w n + c 


(ca + d)w n + c 
w +i= ( ca + d ^ ^ + c 

= w n +2c/(a + d), for n = 0,1,2,..., (1) 

because ca + d = l(a + d) = a — ca. Note that c ^ 0 and 
also a + d ^ 0, because 

(a + d) 2 = (a — d) 2 + 4bc + 4(ad — 6c) 

= 4(ad — 6c) ^ 0. 
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It follows from (1) that w„ —» oo as n —> oo, for all 
Wo € C and hence 

z n = Q+ — -1 a as n -* oo, for all zo € C. 

Remark The sequence w n tends to infinity along a 
straight line, and so z n tends to a along a circle; see the 
solution to Exercise 1.1(c). 

(b) If (a — d) 2 + 46c ^ 0, then / has the two fixed points 
a =Yc { a ~ d + V ( a ~ d ) 2 +46c) 

0=£(a- d -V(*~ d ) a +*bc). 

Putting Wn = h{z n ) = (z n - a)/(z„ - 0), we have 
Z n = (~0w n + a)/(-w n + 1), for n = 0,1,2,... . Hence 
the sequence {z„} is conjugate to the sequence 
«M-1 = h(z n +l) 

= h(f(Zn)) 

= {hofoh~ 1 )(w n ). 

Composing the Mobius transformations h, f and ft -1 (for 
example, by multiplying the corresponding 2x2 
matrices), we obtain 


_ {—a/3 - b + ca/3 + da)w n + (aa + b - ca 2 - da) 
Wn+1 ~ (-a/3 - 6 + c/3 2 + d/3) w n + {aa + b- ca/3 - d/3)' 
Since aa + b = a(ca + d) and a/3 + 6 = /3(c/3 + d), it 
follows that 

_ (<* — /3)(c/3 + d)w n + 0 
n+1 0 + (a-/3)(ca + d) 

-(£$)- . 

/'(«) - 


(ca + d ) 2 ’ 

and so to obtain w n+ i = f'(a)w n , for to = 0,1,2,..., i 
need to check that ad — be = (ca + d)(c /3 + d). Since 

ca + d = 5^a + d+ y/(a — d) 2 + 46c ^ 


and 

c/3+ d = | (a + d - y/(a - d) 2 + 46c) , 
we obtain 


(ca + <f)(c/3 + d) = 1 ((a + d) 2 - (a - d) 2 - 46c) 


as required. 

If |/'(a)| < 1, then w„ -» 0 as n -> oo, for all w 0 € C, 
and hence 



for all 2 0 6 C - {/3}. 

If |/'(a)| > 1, then w n -» oo as n -* 
w 0 S C - {0} and hence 


for all 2o € C - {a}. 


for all 


If |/'(a)| = 1, then |io n | = |tao|, for n = 1,2,..., and so 
the sequence {w n } remains on the circle with centre 0 
and radius |u> 0 |. Hence the sequence {2„} remains on the 
image of this circle under /i -1 , which is a generalized 
circle with a and 0 as inverse points. 


Section 2 

2.1 Using Theorem 2.1, we find that 

2„+1 = 3z n (l - z n ) = -3z* + 32„, n = 0,1,2,..., 
is conjugate to 

Wn+i =wl + d, n = 0,1, 2,..., 
where w n = -3 z n + 3/2, for n = 0,1,2,..., and 
d = -3 x 0 + 3/2 - 9/4 = -3/4. 

2.2 (a) If |c| < \ and |2| < i + y/{ - |c|, then 
|P c (2)|=|2 2 + C | 

< \z\ 2 + |c| (Triangle Inequality) 

< (l + vT^n) +w 

= * + y/k ~ M +1 - M + |c| 

= \ + \J\- W. 

as required. 

(b) If |c| < j and |2| < \ + yj\ - |c|, then, by part (a) 
applied repeatedly, 

l^c n (*)| < 5 + y/\ - w, for n = 0,1,2,..., 
so that Pc(z) -^ooasn-too, and hence z £ E c . 

Therefore, if |c| < |, then we have 

{* : W<i + \A-|c|Jc* c , 

as required. 

(c) If |c| < i, then, by part (b) and Theorem 2.3(a), 

{* : w < \ + - kl | c K c C | 2 : |«| < i + + |c| J . 

Now, if c is close to 0, then \ + y'i ± |c| are both close 
to 1, and so K c is approximately equal to the closed unit 
disc. 

2.3 (a) Since 

f(i) = -i and f(-i) = i, 

the point a = i is periodic with period 2, and belongs to 
the 2-cycle i, —i. Since f'(z) = —1, the multiplier of this 
2-cycle is, by Theorem 2.4, 

f'(i)f'(~i) = (-1) x (-1) = 1. 

Thus this 2-cycle is indifferent, and so i is an indifferent 
periodic point of / with period 2. 

(b) Since 

/(i(-l + v^)) = i(-l- V5) 
and 

/(i(-l- v'S)) = I(-l +V5), 
the point a = j(—1 + V5) is periodic with period 2, and 
belongs to the 2-cycle i(-l + ^5), 1(-1 - -v/5). Since 
f'(z) = 2z, the multiplier of this 2-cycle is, by 
Theorem 2.4, 

/'(|(-1 + a/5)) f'{\(-1 - a/5)) = (-1 + a/5)(-1 - a/5) 
= -4. 

Thus this 2-cycle is repelling, and so |(-1 + a/ 5) is a 
repelling periodic point of / with period 2. 


61 



Section 3 

3.1 


(c) Since 

f(0) = i and /(») = 0, 

the point a = 0 is periodic with period 2, and belongs to 
the 2-cycle 0, i. Since /'(z) = 3z 2 , the multiplier of this 
2-cycle is, by Theorem 2.4, 

/'(0)/'(i) = 0. 

Thus this 2-cycle is super-attracting, and so 0 is a 
super-attracting periodic point of / with period 2. 

(d) Since 
/(e*’/ 13 ) = e 3 **/ 13 , 

and 

/(e 9 "/ 13 )=e 27 * i/13 =e"/ 13 , 
the point a = e n / 13 j s periodic with period 3, and 
belongs to the 3-cycle e* l/13 , e 3lr *^ 13 , e 9lr *^ 13 . Since 
f'(z) = 3z 2 , the multiplier of the 3-cycle is, by 
Theorem 2.4, 

= (3e 2 * i/13 )(3e 6 * i/13 )(3e 18 * </13 ) 

= 2 7e 2<W13 
= 27. 

Thus the 3-cycle is repelling, and so e”/ 13 is a repelling 
periodic point of / with period 3. 

2.4 (a) Since Po(z) = z 2 , a is a periodic point of Pq if 
and only if it satisfies the equation z 2 = z, for some 
positive integer p. Now 



Thus either a = 0, which is a super-attracting fixed point 
of P 0 , or a is a (2 P - l)th root of unity. Since Pq{z) = 2 z, 
the multiplier of a has absolute value 2 9 for some factor q 
of p and a is therefore repelling. Thus the repelling 
periodic points of Pq are the (2 P — l)th roots of unity, for 

P = 1,2,... . 

Now, if T is an arc of the unit circle which subtends an 
angle 9 > 0 at 0, then we can choose p so large that 



and this guarantees that at least one of the (2 P — l)th 
roots of unity (which are evenly spaced around the unit 
circle) lies in T. 

(b) Since Pq(z) = z 2 ", a is a backward iterate of 1 
under Pq if and only if it satisfies the equation 
z 2 " = 1, 

for some positive integer n, that is, if and only if a is a 
2"th root of unity for some positive integer n. 

Now, if T is an arc of the unit circle which subtends an 
angle 0 > 0 at 0, then we can choose n so large that 



and this guarantees that at least one of the 2"th roots of 
unity lies in T. 



The iteration sequence {x n } converges to the indifferent 
fixed point \ of ft/4- 

3.2 (a) First we plot y = x and y = —^ on the same 
diagram, and apply graphical iteration with various 
initial points xq. 



Since y = x lies above y — x/(x + 1), for x > 0, graphical 
iteration shows that if so > 0, then x n —> 0 as n —» oo. 

If x 0 < -1, then xi = X ° , > 0 and so again x„ -* 0 as 
Xo + 1 

n —> oo. If — 1 < xo < 0 and xo A, then graphical 
iteration shows that x no < — 1 for some positive integer 
no and hence x„ —> 0 as n —► oo. 

Since x„ evidently tends to 0 if xo = 0, we deduce that 
x n —* 0 for all initial values xo in R — A. 

(b) The function 



is a Mobius transformation, with a = 1, b = 0, c = 1, 
d = 1, which satisfies (a - d) 2 + 46c = 0. Hence, by 
Exercise 1.4(a), a = 0 is the only fixed point of / and the 
iteration sequence 

z n+ i=f(z n ), n = 0 , 1 , 2 ,..., 
converges to a = 0 for all initial points zo in C. The result 
in part (a) is therefore a special case of Exercise 1.4(a). 
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Section 4 

4.1 (a) Since |P C (0)| = |c| = VS > 2 if c = -1 + 2i, we 
deduce, by Corollary 1 to Theorem 4.1, that c£ M. 

(b) Since P_i(0) = —1 and P-i(—1) = 0, we deduce 
that the terms of the sequence {P" x (0)} are 

0, —1,0, —1,..., and so 

|P-i(0)| <2, fcrn = l,2. 

Hence — 1 6 M. 

(c) For c = — 1 + i, 

Pc(0) = -l + i, 

Pc(0) = (—1 + i) 2 + (—1 + i) = —1 — i, 

Pc( 0) = (-1 - t) 2 + (-1 + i) = -1 + 3i. 

|P 3 (0)| = ’S/TO > 2, 

we deduce, by Corollary 1 to Theorem 4.1, that c ^ M. 

4.2 Suppose that |c| = 2 but c # -2. Then 
|c 2 + c| = |c(c+l)| = |c| |c+l| = 2|c+l|. 

Now |c + 1| is the distance from c to —1, and so 
|c+l|>l. 



Hence |c 2 + c| = |P 2 (0)| > 2 and so, by Corollary 1 to 
Theorem 4.1, c $ M. 

4.3 (a) Since |c| > r c , if |c| > 2, we deduce by 
Theorem 2.2 that if |c| > 2, then the sequence 

|Pc n (c)| = |P" +1 (0)| , n = 0,1,2,..., 

is increasing and so all its terms exceed 2. Hence 

Mu C {c: |c| < 2} = Mi, forn = 2,3. 

(b) Since r c < 2, if |c| < 2, we deduce by Theorem 2.2 
that if |c| < 2 and |P"(0)| > 2, then 

|P c n+1 ( 0 )| = |Pc(P"( 0 ))| > |P"( 0 )| > 2 . 

Thus, for |c| < 2, we deduce that 

|P" +1 (0)| < 2 =» |P"(0)| <2, for n = 1,2,..., 

that is, M n +i C M n , for n = 1,2,..., as required. 

4.4 (a) If c = —1.1 — O.li, then 

|c + 1| = | - 0.1 - 0.1*| = 0.1414... < i. 

Hence, by Theorem 4.4(b), P c has an attracting 2-cycle 
and so, by Theorem 4.3, c € M. 

(b) If c = 0.6i, then |c| 2 = 0.36 and Rec = 0, so that 
(8|c| 2 - §) 2 + 8Rec = (8 x 0.36 - 1.5) 2 
= 1.9044 < 3. 

Hence, by Theorem 4.4(a), P c has an attracting fixed 
point and so, by Theorem 4.3, c € M. 


4.5 (a) If some point a, say, of dlZ lies outside M, then 
because M is closed there is an open disc D with centre a 
which lies entirely outside M. Since TIC M, the open 
disc D does not meet H, and this contradicts the fact 
that a is a boundary point of TL. Hence dH C M. 

(b) (i) From Figure 4.12, it appears that the point 
c = 5 + 5* lies on the cardioid with parametrization 
7 (t) = - \e 2it (i € [-it.jt]). 

To verify this, either note that 

7 (V2)=|i-i(-l) = A + li, 

or check that c satisfies 
(8|c| 2 — |) 2 + 8Rec = 3. 

Thus c lies on the boundary of the periodic region in M 
given by Theorem 4.4(a), and hence c € M by part (a), 
(ii) The point c = — 1 + \i lies on the circle |c + 1| = 
Thus c lies on the boundary of the periodic region in M 
given by Theorem 4.4(b), and hence c € M by part (a). 


4.6 By Theorem 4.5, the function P c has a 
super-attracting 4-cycle if and only if 

P c 4 (0) = 0, but P c (0),P 2 (0) and P c 3 (0) are non-zero. 
Now P c (0) = c, P 2 (0) = c 2 + c, P e 3 (0) = c 4 + 2c 3 + c 2 + c 
and P 4 (0) is of the form 

P 4 (0) = c 8 + 4c 7 H-1- c 2 + c. 

If P„(0) = 0 or P 2 (0) = 0, then P 4 (0) = P 2 (P c 2 (0)) = 0, 
and so P 4 (0) can be factorized as follows: 

P 4 (0) = (c 2 + c)(c 6 + 3c s + ■ ■ • + 1). 

Thus the values of c for which P c has a super-attracting 
4-cycle are the solutions of the equation 
c 6 + 3c s + • • • + 1 = 0, 

and there are at most six such (distinct) solutions. 

Hence, there are at most six values of c for which P c has 
a super-attracting 4-cycle. 

Remark In fact it can be shown that each of the 
polynomial functions c i—> P c n (0) has only simple zeros. 
This makes it possible to count the number of values of c 
for which P c has a super-attracting p-cycle, for each 
positive integer p. In particular, there are exactly six 
values of c for which P c has a super-attracting 4-cycle, as 
indicated in Figure 4.8. 

4.7 For c — — |, P c has the 2-cycle 
(by Problem 2.7(c)) with multiplier 

(^5/O'(«i) = 4a ia2 = 4 (-i + - yfV) 

Since —1 is a primitive square root of unity, we deduce, 
by Theorem 4.6(b), that a period-doubling bifurcation 
occurs at c = —|. 

In Figure 4.8, we see that the point — | lies where a 
periodic region with period 2 and a periodic region with 
period 4 touch, so a period-doubling bifurcation is visible. 


63 




